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SUMMARY 


This investigation is concerned with the interaction 
between a turbulent flow and certain types of structiore re- 
sponding to its excitation. The turbulence is typical of 
those associated with a boundary layer, having a cross- spectral 
density indicative of convection and statistical decay, A 
number of structural models are considered in the investigation. 
Among the one-dimensional models are an unsupported infinite 
beam and a periodically supported infinite beam. The first 
model is used to develops the basic ideas which are then 
applied to the more realistic second model resembling 
the fuselage construction of an aircraft. For the two-dimen- 
sional case a simple membrane is used to illustrate the type 
of formulation applicable to most two-dimensional structures. 
However, a small random variation in the membrane tension is 
included in the analysis since ideally uniform tension never 
exists in practice. Moreover, the mathematical approach used 
in dealing with random membrane tension can be adapted to 
treat other random structural properties in general. Both the 
one-dimensional and two-dimensional structures mentioned above 
are backed by a cavity filled with an initially quiescent 
fluid to simulate the acoustic environment when the structure 
forms one side of a cabin of a sea- or air-craft. 

It is shown that a decaying turbulence can be con- 



structed from superposing infinitely many components, each of 
which is convected as a frozen-pattern at a different velocity. 
This superposition scheme reduces greatly the computation time 
by reducing to one -half the n\imber of integration which must 
be performed on a computer. Furthermore, the scheme provides 
a convenient way in which experimentally measured cross- 
spectral density. of the turbulence pressure fluctuation can 
be incorporated directly in the computation. 

The results of the structure -turbulence interaction 
are presented in terms of the spectral densities of the struc- 
tural response and the perturbation Reynolds stress in the 
fluid at the vicinity of the interface. It is found that 
important spectral peaks of the structural response will not 
appear if decays in the turbulence is neglected in the analysis. 
Thus, the usual Taylor's hypothesis of frozen- pattern turbulence 
is unconservative as far as the assessment of structural 
reliability is concerned. The perturbation Reynolds stress 
is indicative of the change in the skin-friction drag due to 
structural motion. It is shown that, given the statistical 
information of the boxmdary-layer turbulent pressure field, the 
pertvirbation Reynolds stress can be changed Tby varying the 
structural parameters. Therefore, the present study is poten- 
tially useful for designing flight or marine structures to 
minimize the total skin-friction drag. 
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I. INTRODUCTION 

In recent years* there has been considerable Interest 
in the response of panel systems to random pressure fields. 

An aircraft fuselage excited by boundary-layer turbulence or 
by the efflux of a jet engine is a good example. Theobjec-. 
tives of such investigations ZX-93 vary from minimizing struc- 
tural fatigue damage to reducing noise radiation from the 
panel system to either outside or inside of the cabin. 

Recently* another important application has been suggested* 
namely* to design a panel system for the exterior of a vehicle 
such that the total skin-friction drag force over the vehicle 
is a mlnimTim. 

The typical panel system of an aircraft fuselage is a 
multi-span structure which is characterized by close clustering 
of natural frequencies in well-defined frequency bands. The 
usual normal mode formulation does not lead to practical 
results in this case* since it is almost impossible to calcu- 
late the normal modes of structures with a large nximber of 
spans due to close proximity of natiural frequencies in each 
frequency band. However* if all the panels in a system are 
identically constructed* then the structural configuration is 
spatially periodic and the analysis can be greatly simplified. 
Two alternative methods are available to solve the response 
problem of such a spatially periodic structure i the transfer 
matrix technique QlO-16] and the wave-propagation approach 
Cl7->213< The computational simplicity of both methods is 
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achieved by utilizing the fact that the entire system is 
composed of identical subxmits. Although these two methods 
are related [ 222 1 tremsfer matrix technique is more 
suitable for analyzing a periodic structure with finite total 
length whereas the wave-propagation approach is more suitable 
for an infinite periodic structure. 

In general, the spatial periodicity is no longer 
preserved when the excitation is included in the formulation, 
especially when the excitation is a random field. In princi- 
ple, it is possible to express the total response under 
arbitrary excitation in terms of a fundamental solution which 
is the iresponse due to one point-load on the structure. 
Mathematically, this fundamental solution is the Green's 
function and the total response can be represented as a convo- 
lution integral, but the actual calculation can become extreme- 
ly tedious. One type of random excitation which does not 
destroy the spatial periodicity of the system is that which 
is convected as a frozen- pattern at a given velocity 
Known as Taylor's hypothesis, this is an assumption frequently 
made in the analysis of airplane response to atmospheric 
turbulence. Unfortunately, significant decays in the corre- 
lation of pressxire field have been fotond in experimental 
measurements of boundary-layer turbulences C2 3-283. Thus, 
calculations based on frozen- pattern models are just crude 
estimates as far as structural response to boundary- layer 
turbulence is concerned. Recently a scheme has been proposed 
[ 29 - 30 D in which a decaying turbulence field is constructed 
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from frozen-pattern components, thus retaining the computa- 
tional advantage of periodic structure models. This scheme 
(which we shall call the "txirbulence decomposition scheme”) 
will be used in this thesis and applications to structural 
response analyses will be discussed. 

More specifically, the advantages of the turbulence 
decomposition scheme will be discussed in Chapter II by 
comparing to the conventional point-load analysis. In Chapter 
III the basic concepts of the turbulence decomposition scheme 
and its application to the structural response spectriun 
calculation will be developed using a simple model of an 
infinite unsupported beam exposed to a supersonic boundary- 
layer turbulence excitation. The effects of surrounding 
fluid will be taken into consideration. In Chapter IV the 
unsupported beam will be replaced by an infinite beam on 
evenly spaced supports which is a more realistic model of an 
actual aircraft panel system. The solutions will be compared 
with the experimental results. In Chapter V the case of a 
membrane with random structural properties will be considered 
and the effects on the spectral densities of the structural 
responses will be discussed. With potential applications to 
skin-friction drag reduction designs in mind. Chapter VI will 
be devoted to the analysis of the Reynolds stress in the 
botmdary layer. 





II. SPECTRAL ANALYSIS 


2 • 1 Introduction 

Measured cross-spectra of a turbulence field usually 
show some decay in the statistical correlation in addition to 
convection at a characteristic velocity C23-28], Under such 
a random excitation the computation of structural response 
statistics becomes much more tedious than that which would be 
the ease if the turbulence were convected without decay t i.e.« 
convected as a frozen-pattern The conventional method 

of analysis is a point-load approach. As it will be shown in 
the next section, this method requires a numerical double 
integration to compute the cross-spectral density of the 
response of a one -dimensional structure. If the forcing field 
is a convected frozen-pattern field, then an alternative 
formulation will allow the cross-spectral density of the 
structural response to be computed without ntimerical integra- 
tions. However, because of the spatial decay in the measured 
turbulence spectra, the analysis for structural response based 
on the frozen- pattern assumption is just a crude estimate. 

The method to be discussed in this thesis retains the maximum 
computational benefits of the frozen-pattern analysis but at 
the same time the actually measxared spectrum of turbulence 
field can be incorporated in the calculation. In this method 
a decaying txarbulence will be treated as a superposition of 
f rozen- pat tern : components so that the structtiral response can 
be superposed similarly. In the case of a one-dimensional 
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model, this method req.ulr;es oxily a singla integration instead 
of a double integration in the point-load analysis. The method 
will be called the "turbizlence decomposition scheme" in the 
sequel. 

To extend the turbiilence decomposition scheme to a two- 
dimensional problem, the actual turbulence field is divided 
into strips which are running parallel to the x-eixis, the 
direction of the mean flow. Each strip can then be decomposed 
into infinitely many frozen- pattern components as in the one- 
dimensional case. The total structxirsLL response is obtained 
by summing up the responses due to all the loading strips. 

In essense the point-load scheme is used in the across-flow 
direction whereas the frozen-pattern scheme is used in the 
along-flow direction. 

In the following sections the one-dimensional case will 
be discussed first and then the two-dimensional structure. 
Furthermore, a general treatment of random structural properties 
will be included in the formulation of the two-dimensional 
problem. 

2.2 The Statistical Properties of the Excitation 

Measured frequency cross-spectra for pressure varia- 
tions in a turbulent boundary layer with respect to a fixed 
frame of reference have the general form of 

* ip(0,0,&i))il;i(?)i|f2(T)exp(-i(i)?/Uc) 


( 2 . 1 ) 
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where U is a characteristic velocity* and t|r and iju are non- 

C X » 

negative definite even functions of % and , respectively. 

These two functions have an absolute maximum equal to one at 
the origin and they approach to zero at large absolute values 
of the argument. The general form* Eq. (2,1)* is sometimes 
attributed to Corcos [ 23 J, A number of researchers have 
reported curve-fitted results for l|p(0*0*u3)* and 

For representative works we cite the papers by Bull 
Willmarth and Wooldridge [|25ll* 2 Uid Maestrello* et al [[26-283, 
Implicit in Eq. (2.1) is that a real turbulence is not a frozen 
one. We note that the above pressure spectrum reduces to that 
of a plane wave field if ~ and it would correspond to 

a frozen-pattern turbulence if = 1. 

For additional physical insight and later use in sample 
calculations* two measured spectra are given belowt 


For a subsonic boundary layer C26*273» 


$ (0*0*co) 
P 


if 

2 U« 


3 

z 

n=l 


An 



( 2 . 2 ) 


\|jj^(e) » exp ( 



(2.3) 


Aj = 0.240 = 0.470 

Ag = 1.08 K2 = 3.0 

kj = 1.80 = 14.0 

= 2 
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-3 -3 

0 * - (1,24 X 10 ) (I^ysi ) + 1.15 X 10 sec 


0.8 U„ 


where Um is the free-stream velocity, denotes the speed of 
soxmd in the field, S* represents the houndary-layer displace- 
ment thiclcness, and 0 is the average eddy lifetime. 

For a supersonic boundary layer E28Di 


— 1 s ^ 

*«(0,0,00) = - ~ S A_ 
♦j(5) =. exp ( - ) 

♦2(7) - exp ( 

- ^ Knl"i 

= 4.4 X 10“^ 

-2 

= 5.78 X 10 

-2 

Ag = 7.5 X 10 

Kg = 2.43 X 10"^ 

* -9.3 X 10“^ 

= 1.12 

A4 = -2,5 X 10“^ 

K4 = 11.57 

= 3 

o<2 = 0.26 


(2.5) 

( 2 . 6 ) 

(2.7) 


Uc = 0.75 Uoo 


Where S is the boundary-layer thickness which is defined as 
the distance from the boimdary at which the average value of 
the turbulence velocity reaches 0,99 Uo,. 
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2.3 Foin-t>Load Analysis 

Consider a unit concentrated load at x » ^ on a one- 
dimensional structure where x is measxired along the structwe. 
parallel to the flow direction. We obtain a frequency response 
function “HCx, I ,u>) by solving the equation 

.C C H(x,§,to) e^“^ } * S(x-6) (2.8) 

where JLC } represents a linear differential operator in x 
and t, and S( ) is a Dirac delta function. Then, the cross- 
spectral density of the structural response iXgfCo) is 

calculated from 

■ff- 

= JJ H(Xj^,§^,cu) H*(x2,?g,w) 

00 

d^idSg (2.9) 

where $p(Sj**52*^J cross-spectrvun of the pressure field, 

i is the length of the structure under consideration, and an 
asterisk denotes the complex conjugate. It is assumed that 
the loading is a spatially homogeneous random process t there- 
fore, the cross-spectrum $p depends on Although this 

method is quite straightforward, long computer time is required 
to carry out the double integration in Eq. (2.9). 

2.4 Frozen- Pattern Analysis 

If the presswe is truly of a frozen type and is 
convected at a constant velocity in the positive x-direc- 
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tiont then it is a r^dom fimction of x - U^t'. Such a random 
function can he expressed as a Fourier-Stieltjes integral as 
follows t 

p(x - Ugt) » I e^^“^ ” dP(]c) (2.10) 

itliftre < the frequency o) and the wave-nxunber Ic are related to the 
convection speed as cc/h It is known from the random 

process theory that 

ECdPCki) dP*ac2)5 “ Sp(ki) SOcj-kgX diidkz (2.11) 

where E£ } represents the ensemble average, and Sp(k) is the 
wave -number spectrum in a coordinate frame moving at the 
velocity I7 q (referred to as the moving frame in the sequel). 

The cross-correlation function ECp(xj^-Ugtj^ ) p(x 2 -UQt 2 )} 
of the pressure, referred to the fixed frame, can be calculated 
simply by use of Eqs. (2.10) and (2,11), This function, 
denoted by Rp, depends only on £ - U^t: where § = x^ - Xg and 
t: = tj - t 2 » and it is related to the moving-frame wave-number 
spectrum Sp(k) as follows i 

r ik(Uc-c - 5 ) 

Rp(§ - U^jT) = J e Sp(k) dk (2.12) 

— oo 

If a Riemann-Pourier transform is taken of Eq. (2.12) we 
obtain the fixed-frame frequency cross-spectrum of pi 



10 




^ a¥ j 




1 „ ,c0x -i(^/Uc)5 

ss s (•^^1 e 

|Ucl 


(2.13) 


Equation (2,13) shows that the fixed-frame frequency cross- 
spectrum of a frozen- pat tern txarbulence has a special form 
where 5 appears only in the imaginary exponent. This equation 
also provides a simple formula to convert from S (k) to 
^_(5»0)). Conversely, to convert from $_($,W) to S (fc)» 


Sp(k) =|uj$p(0,ku^) 


(2.14) 


Evaluated at ^ = 0 the cross-spectrum $p($,CO) reduces, of course, 
to the usual spectrum. We emphasize that Eqs,. (2.13) and (2.14) 
are valid only if the turbulence is strictly of a frozen- 
pattern, and is convected at speed U^, 

Equation (2,10) suggests that the structural response 
to a frozen- pattern turbulence can be constructed from a 
fundamental solution where the excitation is just a convected 
sinusoidal pattern of unit amplitude. Thus, let H(x,k) exp(ia)t) 
be the steady-state solution for 


-CCH(x,k) exp(iwt)} = expCi(£*3t - kx)] (2.15) 


Of course, this solution must satisfy all the necessary 
boundary conditions. Then the solution to 

«C-Cw(x,t)} = p(x - Ugt) 


( 2 . 16 ) 
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after reaching stochastic stationality, may he expressed as 

oo 

w(x»t) = J H(x»k) exp(ioot) dP(k) 

— ao 

OO 

= I H(x,k) expdUjjkt) dP(k) (2,1?) 

— OO 


It follows that the cross-correlation fimction of the structural 
response is 


ECw(Xj^,tj^) wCXgftg)} 

(f iUf» (^1 1.« — kpto ) 

— OO 

Sp(lCj^) SO4 - kg) dkjdkg 

• "t ) 

= J H(Xj,,k) H*(x2pk) 2 Sp(k) dk (2.18) 


As expected, this correlation function is dependent only on 
tf - tg. If it is desired to calculate this correlation func- 
tion in the frequency domain, we may substitute Eq. (2,l4) 

Into Eq. (2.18) and change U^k to (Oi 

E{w(Xj^,tj^) w(Xg,tg)3 

( , iw(t, - tp) 

= J H(xj^,a>/Ug) H*(xg.cu/u^) e ^ ^ $p(0,(jO)dio (2.19) 


In terms of the input and output spectra the relations are 
extremely simple and illtuninative i they aret 


12 


in the wave -number domain i 

S^(xj,X2»k) = H(Xj^,k) H*(x2fk) Sp(k) (2.20) 

in the frequency domain i 

$^(Xj, ,X2»oi) = H(Xj^.o)/Uq) H*(x2.W/Uq) $p(o,o)) (2.21) 

When » X2 these formulas reduce to those for the usual 
spectra » and they have the same form as the well-known result 
for a single degree of freedom system in the random vibration 
theory. 

It is appropriate to call the H function in Eqs. 

(2.20) and (2.21) the wave-number response function for 
convected frozen load to distinguish it from the K function 
in Eq. (2.9) which is the frequency response function for 
point load. 

The advantage of the frozen-pattern assumption is 
clear. To obtain the cross-spectrtim of the response, no 
integration is required under this assxmption while a double 
integration is needed in the point-load analysis, Eq. (2.9). 
Unfortunately, significant decays in the correlation have 
been found in experimental measurements of boundary-layer 
turbulences. Thus, correlations based on frozen-pattern 
models are just crude estimates, at least for the calculation 
of structural responses to the boundary- layer turbulence. 
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2,5 Turbulence Decomposition Scheme - One-Dimensional Case 


In the one -dimensional case measured frequency cross- 
spectra have the general, form of 

*p(5i‘*>) “ ip(0,w) exp(-iw5/UQ) (2.22) 

which is reduced from Eq. (2,1). Implicit in Eq. (2.22) is 
that a real turbulence is not a frozen-pattern one. 

To obtain a theoretical spectrum consistent with Eq. 
(2.22) the following representation of a general turbulence 
pressure is proposed t 

OO 

p(x,t) = ( p(x - ut) dG(u) (2.23) 


Eq, (2,23) implies that p(x,t) is a superposition of infinitely 
many frozen- pattern components, each having a real random 
amplitude dG(u) and a convection velocity u. Such velocities 
cam assume either positive or negative values. Of course, 
each frozen-pattern component can, again, be decomposed into 
sinusoids. Thus 



i(ukt - kx) 
e 


dP(k,u) dG(u) 


(2.24) 


and its fixed-frame correlation is 

E{p(xj^,tj^) p(xg,t2)} 
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= Jill expCKuj^kltj^ - U2k2t2) - - ^2X2)3 

— 00 

ECdP(kj^,Uj^) dF*(k2,U2) dG(u^) dG{u2>} (2.25) 

In order that this correlation function may depend only on 
^ = Xj^ - X2 and “C = tj^ ” ■^2* which we shall assume to be true, 
the ensemble average under the integral sign in Eq. (2,25) 
must have the form 

ECdP(kj^,Uj^) dP*(k2,U2) dG(u^) dG(u2)3 

= Sp(kj^,Uj^) S(kj^ - ^ 2 ) S(uj^ - U2) dkj^dk2duj^dU2 (2.26) 
Substitution of Eq. (2,26) into Eq. (2.25) results in 

00 

Rp(?,r) = II e^^^^ " Sp(k,u) dk du (2.27) 

— 00 

We now apply a Fourier transformation to obtain the fixed- 
frame frequency spectrum 

T 

- Sp(f.u) du (8.28) 

— QO 

Clearly Eq, (2,28) is a generalization of Eq. (2,13), 

To compare Eqs. (2.28) and (2.22), the latter is 
Fourier-transformed to yield 



where 


^(v) * 2^ j (2.30) 

— oo 

Therefore » 

CO 

fp(5,W) = $p(0,a3) I - -g^) dot (2.31) 

— oo 

Letting c^ = oyu, we obtain 


$p($,C0) 


$p(0,cc) 


00 

1 


CO 

2 


lu 


?(“ - f-) 


dn^ (2.32) 


Then equating $p and $p we find a f ormtila to compute 
Sp(u>/u,u) as follows t 


' ItI - -5;) v°’“> 

Oo 

= Is^l $p(0,w) j 4f(5) expCi?(^ - ^)1 (2.33) 


The frequency cross-spectrum for the structural 
response can be obtained by a similar superposition* Thus 
by a generalization of Eq. (2.21), 




■Jii 


H(Xj^,^) H*(x,,^) S„(^.u) du (2.3^) 


"2»u" P'u 


Or, letting k = 



16 


OO 

= J |5jH(Xj^,1e) H*(x2fk) SpCk,^) (2.35) 


NoWf since 


Sp(k,^) = |k| $ (0,oj) f (k -^) 

c 

we obtain a very simple result 

Oo 

^^(Xj^.Xgiw) = ip(0,U)) j H(x^,k) 

— Op 

H*(x2.k) ¥(k - ^) dk (2.36) 


As a check we note that when the turbulence is frozen-pattern, 
tj;(S) = 1 and 

?(v) = S(v) (2,37) 

then Eq,. (2,36) reduces to the same form as Eq.(2.21), 

In the case when the structural acceleration response is 
the main concern of the problem, Eq, (2,36) is changed to 

OO 

$^(Xj^,X2»o^) = a? ip(0,U3) j H(Xj^,k) 

— OO 

H*(x2,k) Y(k - ^) dk (2.38) 

2.6 Turbulence Decomposition Scheme - Two-Dimensional Case 

Generalization of the results obtained in the last 
section to a two-dimensional, pressiure field is straightforward. 
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Assuming that such a pressure field can be decomposed Into 
frozen-pattern components » we have 

p(x - ut,y) dG(u) (2,39) 

where both p and G are reuidom functions and $ is a frozen- 
pattern component of p. Each component pressure ^ can again 
be constructed from frozen-pattern sinusoids. Thus 


p(x.y.i;) 


I 


p(x,y,t) 



i(ukt - kx) 
e 


dP(k,Uty) dG(u) 


(2.40) 


In order that the cross-correlation of the random pressure 
ECp(x^fyj^ttj^) P(x2»y2»t2>3 1 b dependent only on S = Xj^ - x^, 
7 “ yi ” y 2.* = tj^ - t2» the cross-correlation of dFdG 

must have the form 

E£dP(kj^*u^,yj) dP*(k2»U2,y2) dG(Uj^) dG(u2)} 

= Sp(kj,Uj,T^) S(k^ - k2) S(Uj^ - Ug) 

dkj^ dk2 dUj^ dU2 (2,4l) 

Then, it can be shown that the cross-spectrum of p has the 
form 

7 1 -iwS/u 

$p(5,'7f<^) =J juj ® (2.42) 

— oo 

Equating Eq. (2,42) with the general form of the measured 
turbulent x^essure spectra, Eq. (2.1), results in 



IS 


Sj,(g,u.^) - - f.) 5j,(0.0.w) teO[) (2.43) 

c 

where 'Pj^ is the Fourier transform of 1 };^^ i i,e, , 

oe 

’Pl(v) = ^ j (2,44) 

— 00 

It is clear from Eq, (2.40) that, in this case, the 
structural response can also be constructed from a fundamental 
solution which is obtained by letting the excitation be some 
suitable frozen-pattern sinusoid. Denote this fxmdamental 
solution by H(x,k,y,y* ,o)) and let H exp(iwt) be the steady- 
state solution for 

i(jot i(wt - lex) 

XCH(x,k,y,y’,w) e 3 = e S(y - y') (2,45) 

where JC, represents a linear operator in x, y and t. The 
physical meaning of H is self-explanatoryj it is the complex 
amplitude of the steady-state structural response at (x,y) 
due to a strip of excitation along y = y* which is convecting 
in the x-direction. Then, the solution to 

»CCw(x,y,t)} = p(x,y,t) (2,46) 

after reaching the stochastic stationality, may be expressed 
as 

b ^ 

w(x,y,t) = j 

0 —00 

dP(k,^,y) dG(-^) 
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b 

JdyJj 


H(x,k,y,y* ,ku) e^^^dF(k,u,y' ) dG(u) (2.^7) 


0 

where h is the width of the structure in the y-direction. The 
cross-correlation of the response is* therefore, obtained as 
follows I 


E£w(xj^,yj^,tj^) w(x2.y2»'t2)} 
b 

II 1 1 

H(Xj^,k^,yj^,y^,k^Uj^) 


b ^ 

= j\ \[ 


, - i(k. u. t^ - k_u_t_) 

H*(x 2 ,k 2 ,y 2 .y|.k 2 ^ 2 ^ e 1 1 1 222 


Sp(kj^,Uj^,y£ - y-|) S(kj^ - k 2 ) S(uj^ - Ug) 


dk^^dkgdUj^dUg 

b 


= j| dy£dyjJ H(Xj^,k,yj^,y£,ku) H*(x2fk,y2,y^,ku) 


ikui: 

e Sp(k,u,y^ " 


(2.48) 


Substituting Eq. (2.43) into Eq. (2,48) and applying a Fourier 
transformation, we obtain the cross-spectrum of the structural 




= $p(0,0,w) 



0 


dy^dy^ 


H(Xj^,k,yj^,y£,u>) 


response 
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H*(x2,k,y2,y^,co) ?j^(k - - yp dk (2.49) 

c 

The reason for not using the same decomposition scheme in the 
y-direction as that for the x-direction is the absence of true 
trend of convection which can be seen in Eq. (2,1), Therefore, 
the correlation length of the random forcing field in the. y- 
direction seldom extends beyond one panel and the structural 
response can be computed quite accurately using a structxiral 
model consisting of just one row of panels running in the x- 
direction. For such a model the usual Levy series represents-, 
tion of the y-direction response is adequate and the double 
integration on y£ and 73 Eq, -(2,49) .can be carried out 
without difficulty. 

Thus far, the structural properties are assumed to be 
deterministic; therefore, the frequency response function H 
is a deterministic fvinction. Sometimes it may be of interest 
to include the effects of random variability of the structural 
properties in the calculation, then H becomes a random process. 
In this case the cross-correlation of the structural response 
becomes 


E£w(xj,yj^,tj^) w(x2,y2it2)} 


dy^dy^ 


E£H(x^,k^,y^,y^,k^u^) 

H* (x2 , k£ , y2 1 yj » ^ dF(k^ ,Uj^ »y^ ) dF*(k2,U2,y2) 


dG(Ui) dG(U2)3 


(2.50) 
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If the structiiral properties and the turbulence properties are 
independent of each other » which is a reasonable assvunption, 
the ensemble average inside the integral of £q. (2.50) is 
separable. Again, using £q. (2.4l) we have 


E£w(xj,yj^,t^) wCxg.yg.tg)} 



E£H(Xj^,k,yj^,y^,ku) H*(x2ik,y2,yj|,ku)} 


Sp(k,u,y^ - y^) 


ikut: 

e 


dk du 


(2.51) 


The cross-spectrum of the structural response is obtained by 
a Fourier transformation of Eq. (2.51), resulting in 


$^(x^,X2,yi,y2*‘*)) 




b 

,0,co) fj dy£ 


dy^ 


J 


ECH(x^,k,yj^,y£,a)) 


H*(x2,k,y2,y^,uj)3 ij/j^(k - ~) (2.52) 

c 

2.7 Conclusion 


The turbulence decomposition scheme has been discussed 
and compared with the conventional point-load approach. It 
has been shown that a decaying and convecting turbulence 
pressure field can be constructed from frozen- pattern components, 
each having a different convection velocity, and that this new 
scheme simplifies greatly the analyses of random properties of 
structural responses under the excitation of boundary- layer 
turbulence pressure. In paticular, only a single integration 
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is required to compute the cross-spectrum of the struct\a:al 
response instead of a double integration in the point-load 
approach. The integration involves a wave-number response 
function which is the response of the structure to a unit 
convected sinusoid. In the following chapters applications 
of this scheme will be discussed for some structural configu- 
rations and flow field problems. 
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III. UNSUPPORTED INFINITE BEAM 
3 • 1 Introduction 

A simple model which » nevertheless, retains most 
import^t fea,tures of aircraft panels in a boundary-layer 
environment is the infinite beam shown in Pig, 3.1. The 
beam is backed on the lower side by a space of depth d which 
is filled with an initially quiescent fluid of density an.d 
sound speed a. 2 » On the upper side the beam is exposed to the 
excitation of a supersonic boundary-layer turbulent pressure 
p. The fluid on the upper side of the beam which carries 
the turbulence has a free-stream velocity U*,, density and 
sound speed a^^. An actual panel system of .. an aircraft is 
reinforced by stringers and frames so that this unsupported 
beam is not a good representation at low frequency range. 
However, at high frequencies, the turbulence eddy size is 
much smaller than individual panels. Then the effect of the 
constraints at the stringers and frames becomes negligible. 

In any case the infinite beam model is an ideal one which can 
be used to show the utility of the ttirbulence decomposition 
method as well as the structure -fluid Interaction without the 
burden uf mathematical complexities, 

3.2 Wave-Number Response Function of the Infinite Unsupported 
Beam 

As the beam responds to the excitation pressvire p(x,t) 
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its motion will generate additional pressiires in the fluid 
media on the upper and lower sides. Denoting these generated 
pressTires by pj^ and P 2 # respectively, the governing equation 
of the beam motion is given by 


El 




2 

:afw 




<3.1) 


Where £ denotes the Young's modulus, I is the moment of 
inertia, m is the mass per unit length of the beam. 

For the purpose of determining the wave-number 
response function, H(x,k), the turbulent pressure p should 
be replaced by expCi(‘ 2 >t - kx)D and the structural response w 
equated with H(x,k) exp(itot) = A(k) exp(-ikx) exp(ioit). 
Furthermore, we shall make the usual approximation that Pj^ 
can be calculated without regard to the presence of the 


turbulence. Then Pj^ is governed by the equation 

Pi = 0 


(3.2) 


and subject to the conditions that pj^ can propagate only in 
the positive z-domain, and that 




^>2 


i^iot - kx) 


(3.3) 


The solution for Pj^, when evaluated at z = 0, is known to be 
C31] 

,2 


\ \ p a k(o/k - Uoo) 

< Pi>,=0 - 


A(k) e 


i(tot - kx) 


(3.4) 
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Some comments about Eq, (3.^) are in order i (1) o)/k is the 

speed at which the structvural motion A(k) expCi(uit - kx)3 is 

propagated along the beam, (2) A structiiral motion generates 

no pressure in the adjacent fluid medium if it is propagated 

at the same velocity as that of the fluid medium (the case of 

w/k = Uoo)t (3) Theoretically, the generated pressure attains 

an infinite amplitude when the propagation velocity of the 

structural motion relative to the medixim is equal to the speed 

of sound (the. case of - Uool = the shock-wave effect), 

(4) When this relative velocity is less than the speed of 

soundt i.e,, |^/k - Uo.| <a^, the generated pressure should 

provide additional inertia for the structiu*al motion (the 

apparent mass effect); therefore, a negative imaginary value 

2 2 4 

should be given to the square-root C(aj/k - Uoo) - a^^ 3 in the 
calculation. 

The pressure generated on the lower side of the beam 
is governed by the equation 


T:r - ^2 (7-2 + P2 

3t* Sx 3z 

and subject to the conditions 

ap2 


» 0 


dZ 


= 0 


at z = -d 


(3.5) 


(3.6) 


and 


It 

az 


A(k) 


at z = 0 


(3.7) 


The solution for P2» when evaluated at z = 0, 


is given by 
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/ X _ .2 cot(Yd) _ . i(u)t - Tsx) 

(S2)j=0 ■ fsf“ X *<''> ' 


( 3 . 6 ) 


where 




(3.9) 


2 

For small d (shallow cavity) and Y > 0, this pressitre provides 

additional stiffness on the structural motion hut for certain 

ranges of d value it can change to w added mass. When 
2 

Y < 0, V becomes imaginary in which case 


cot (yd) coth(yd( 

^ ‘ lr| 


( 3 . 10 ) 


Again, the j >2 term has the effect of an added. mass on the 
structural response. 

Eqs. (3.^) and (3.8) can now be substituted into Eq. 
(3.1) to find A(k) and. therefore. K(x.k). recalling that p 
must be replaced by expCi(wt - kx)3 and w by A(k)expCi(«Jt-kx)3. 
The result may be expressed as 


H(x.k) 


A(k) e 
4 


-Ikx 


2 k((o/k - U,o)‘ 

“ tEIk - mto + lfiai 5 5^ 


P2" 


2 oot(rd) ^-l -ilQC 
^ j e 


(3.11) 


This wave-number response fxmction H can now be substituted 
into Eq. (2.36) to obtain the cross-spectral density of the 
structural displacement response. 
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3*3 Numerical Example 

Ntmerical computations have been carried out for the 
frequency spectrum (i.e.t when x^ ~ ^ in £q. (2.36) ) of the 
structured, response using the following physical datai 
properties of the beami 

il 9 

El .(bending rigidity) « 3»935 x 10- N-m^ . 
m (isass per unit length) » 9*7^6 Kg/m 
properties of the surrounding fluid media i 

= P2 “ f density) = 0,11015 Kg/m^ 
a^ s a2 ~ a (speed of sound) = 26I.6 m/sec 
Ito (free-stream velocity on upper side of beam) 

= 575.6 m/sec 

d (cavity depth) = 0,1178 m 
properties of the supersonic bovindary- layer turbxilence 
pressure C28]i 

$_(0,w) = spectral density * "r 2 A_ 

i|f(5) * decay factor = exp (- -J^) 

U^\ (characteristic convection velocity of the 
turbulence). ’= 0,75 Uoo 

y(k - =s 5-^ 5— 

Uc TT^SCC'^S)"^ + (lE - o/Uq)^] 

S (boundary- layer thickness) = 0,279 m 
experimentally determined constants 
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II 



i.'2 

Ai » 4,4 X 10 

% 

» 5.78 X 10 “^ 

Ag = 7.5 X 10"^ 

K2 

= 2.43 X 10 “^ 

A3 » -9.3 X 10“^ 

K3 

= 1,12 

-2 

Ajij, * -2,5 X 10 

K4 

» 11.57 


Pig* 3*2 shows the computed displacement frequency 
spectrum of the structural, response imder the assumption of 
a truly frozen->pattern turbulence (i|;(^) - 1). There appears 
cnly one peak aroxind 4?5 Hz and the value of spectrum decreases 
rapidly as the frequency Increases. The existence cf the peak 
can be explained as follows t In the absence of surrounding 
fluids, the wave number of the free structural motion would 

2 i 

be (mto/EI) . The wave number of the frozen- pattern turbulence, 
k», is related to the circular frequency O) by k = A 

-resonanca. (called coincidence) occurs when these two wave 
numbers are equal In the present case this coincidence 

frequency is found at 466 Hz. Therefore, the peak in Fig. 3,2 
results from the coincidence of the wave numbers of the frozen- 
pattern turbulence and the free structural motion. The small 
difference between the peak frequency in Fig, 3*2 and the 
estimated value above comes from the effect of the svtrrounding 
fluids on the structural response. 

Figk 3*3 shows the spectrum of the structvu^al displace- 
ment response when the measiired spectrum of the turbulence 
pressure field is used in the computation. This spectral 
density has many peaks in contrast with only one peak in the 
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frozen-pattern case. Since infinitely many frozen- pattern 
components are present in the turbulence, the turbulence 
spectrum has also infinitely many wave-number components 
at each frequency. Thus, the coincidence (if we may use 
this terminology in this case) occurs at each frequency. 
However, the response magnitude may be small if a frozen- 
pattern component which generates a coincidence wave in the 
structure has little contribution to the turbulence spectrtun 
at that frequency. Therefore, even the coincidence is present 
at each frequency, its magnitude vsiries as the frequency 
changes and some of them appear as pealcs in the response 
spectrum, 

A comparison of Pigs, 3»2 and 3t3 shows that the 
frozen-pattern assumption is xmconservative which leads to 
lower estimates for the structural response and the radiated 
noise level 

3,4 Conclusion 

The theory developed in Chapter II has been applied 
to the simple example of an unsupported beam exposed to 
boundary- layer excitations. The effect of a cavity and the 
effect of the free-stream velocity are included in the 
analysis. The spectral density of the structural displace- 
ment response was calculated using the measured turbulence 
presstire spectrum of a supersonic boundary layer auid the 
result was compared with the solution obtained for the ideal 
frozen- pattern txjrbulence. This comparison has shown that 
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the frosen-pattern assmnptlon leads to considerably lower 
estimates for the structural response. Thus* we may conclude 
that the frozen-pattern assmnption should not be used in the 
computation of the structiaral response spectrum under the 
excitation of a boundary-layer turbulence pressiire. 
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IV. PERIODICALLY SUPPORTED INFINITE BEAM 
4 . 1 Intr oduc ti on 

In this chapter the ttirbulence decomposition scheme 
will he applied to the analysis of an infinite beam supported 
at uniformly spaced intervals by elastic springs. The elastic 
supports are simplified versions of reinforcing stringers of 
an aircraft fuselage. Although an aircraft fuselage is a very 
complicated multi-panel system its dynamic behavior is similar 
to that of the periodically supported beam described above. 

The one-dimensional beam problem, however, is more suitable 
for fundamental studies since basic concepts can be developed 
without the burden of mathematical details. Thus, the analysis 
of the present chapter will be restricted again to one spatial 
coordinate. 

At the first sight the problem of a periodic beam may 
not appear more difficult than that of any other structure if 
one accepts the linearity assumption and uses a normal-mode 
formulation. In practice, however, the normal mode of a 
periodic beam of many spans cannot be calculated accurately 
due to close clustering of natural frequencies in frequency 
bands. The futility of the normal mode approach in dealing 
with a large number of spans has led to two alternatives i the 
wave propagation approach (space -harmonic analysis) Cl 7-21 3 
and the transfer matrix approach ClO-163* The two alternative 
methods are closely related, however. The so-called free wave 
propagation constants in the first method are the natural 
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logarithms of the eigenvalues of the basic transfar matrix in 
the second method £22 J, The computational simplicity in both 
methods is obtained by utilizing the fact that the entire 
system is composed of identicsil sub-units in the formulation. 
The fundamental solution required for the construction of the 
total structural response is one corresponding to the excita- 
tion of a frozen- pattern Sinusoid. To obtain this fundamental 
solution the formulation will follow Mead's wave propagation 
method* but will take into account the effect of free-stream 
velocity on the same side of the turbulence excitation and the 
effect of a cavity on the opposite side of the excitation. As 
a numerical example* the spectral density of the structural 
response will be computed and the results will be compared with 
experimental measurements . 

4.2 Wave-Number Response Function of the Infinite Periodic 
Beam 

A sketch of the structural model is shown in Fig. 4.1* 
surrounded by the same acoustic envircnment as that assxmed in 
Chapter III. 

The governing equation of the beam motion not directly 
over an elastic support is given by 

a^w a^w 

“TJ “ P ■ P2>z=-d 

OX OX 

where D denotes the bending rigidity and m is the mass per 
unit length of the beam. The additional pressure fields 
denoted by pj^ and p£ are generated in the fluid media on the 



p ( turbulence pressure ) 



Fig, 4,1 An infinite periodic beam under the excitation 
of boundary-layer turbulence 
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upper and lower sides* respectively* due to the beam motion. 
For the purpose of determining the wave-nomber- response 
function* H(x*k)* the turbulent pressure p should be replaced 
by expCi(j:ut - kx)3 and the structural response w by H(x*k) 
exp(iurt). Thus* Eq.. (4.1) becomes 

(D-^- mo?H) (p - pg) (4.2) 

&x^ Z“d 


The forcing function expCi(tot - kx)3 gives every span the same 
excitation but with a phase-lag jJq * ki from one span to the 
next. In this sense jXq may be considered as the imposed phase 
lag of the excitation. To satisfy the spatial periodicity in 
the structural response Mead suggested the following series 
form Cl83« 


H(x,k) = Z An exp(- 

n<=-oe » 


(4.3) 


where 


+ 2nTT = kJ + 2nTT 


(4.4) 


Without the elastic supports the wave-nvunber response function 
would be just the one term associated with the forcing phase- 
la^yt^. The elastic supports give rise to multiple reflec-* 
tions* thereby admitting other values. 

For the induced additional pressure p^* we apply the 
same assumption used in the preceding chapter. Then p^ is 
governed by the equation 
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(^. 5 ) 


and subject to the conditions that can propagate only in 
the region z > d, and that 


r. 3 »2 iu>t 


(4.6) 


The solution for p^^ corresponding to each component in the 
wave-number response function exp(-iyijj^x/t) is known. The 
total Pj^ can then be obtained by superposition. This pressure, 
when evaluated at z = d, is given by C3l3 


I - • (u^ - 

^l|z=d n^-oo'^ - Vcc)^- a^^]* 

i((ot - j^j^x/Ji) 


(4.7) 


where 



(4.8) 


The same comments which were made previously relative to Eq, 
(3.4) apply in this case to each term in Eq. (4,7). Each 
component now has a different propagation speed u^ instead 
of the one propagation speed o)/k in Eq, (3.4), For emphasis 
and clarity these comments are restated as follows i (1) is 
the speed at which the component structural motion Aj^ expCi(u)t 
- JJ^x/JL)3 is propagated along the beam, (2) A component 
structural motion generates no pressiire in the adjacent fluid 
mediiim if it is propagated at the same velocity as that of the 



fluid medium (the ease of Uj^ = U«|)« (3) Theoretically, the 

generated pressure attains an infinite amplitude when the 

propagation velocity of the structural motion relative to 

the medium is equal to the speed of soxuid (the case of |u^ - 

Uool = the shock wave effect), (4) When this relative 

velocity is less than the speed of soundi i,e,, jUj^ - Uoo|<aj^, 

the generated pressture should provide additional inertia for 

the structural motion (the apparent mass effect)* therefore, 

a negative imaginary value should be given to the square-root 
2 -2 i 

C(Uj^ - Uoo) - 3 in "the calculation. 

The induced pressiire P 2 in the fluid medium 0 z < d 
is governed by 


2 2 

^ ?2 - 2, a .a > _ 

— - ^2 ( — 5 + — -P; 
bir bx bz ' 


= -o 


and subject to the conditions 


(4.9) 


Jip_ 

"52 ® ® at z = 0 

^P2 ^ 2 i<»>t 

* 5 ^ = ^ ® at z = d 


(4.10) 

(4.11) 


The solution for P 2 , when evaluated at z = d, is given by 



P2“^ 


eo 

z 

n=-<» 


cot(rT,d) ijyt -jd^x/1) 


*n 


'n 


(4.12) 


where 
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2 

For a small d and a positive » the n-th component of this 

transmitted pressiure P 2 gives rise to additional stiffness 

on the structural motion and for certain ranges of d value 

2 

it can become an added mass to the system. When X is 

n 

negative} becomes imaginary in which case 


cot ( Tjjd ) eoth)yy^d| 

° ■ Knl 


(4.14) 


and the component always contributes to the system inertia 
regardless of the value of d. 

Equations (4,3) » (4,7) and (4,12) can now be substi- 
tuted into Eq, (4,2) to obtain 

^2 ^Aj^cp(n) exp(-ijinX/i) = exp(-ij<Q3^) (4.15) 


where 


if 2 

cp(n) = D (-J-) - mo> + ifj^a^ 


(Uj^ — Uflo) 






2 cot(ynd) 


TT 

U„) ~ a-i J 


(4.16) 


'n 


To determine the amplitude we follow Mead*s 
procedure Ci93 and calculate the virtual work done by the 
external forces acting on the structvire and by the internal 
forces in the structure through a virtual displacement 



^-i(cot - 



Excluding the elastic supports, the virtual work done within 
one span of the heaia is 

sw„ - <p(n) J d* 

“0 

0 

The elastic supports are characterized by a translational 
spring constant K^, a translational inertia M, a torsional 
spring constant K^, and a torsional inertia I. Thus the 
virtual work contributed by each elastic spring is 


SW^ = 


SW^ + S\ 


(K. - Mw ) SA. Z A- + (kr. 
^ Jns-oo ^ 

^ 


Io>^) 


(^. 18 ) 


Since the structural motion is spatially periodic the virtual 
work done throughout the entire structure is proportional to 
that of a periodic unit. Therefore, the principle cf virtual 
work can be stated for a periodic unit as follows: 


SWtj + SW^. + SWp = 0 


(^. 19 ) 


which leads to the simultaneous algebraic equations: 

2 2 

, ^ K+ - Mo) oo Kr - 1(0 oo 

AjCp(j) + ^ ^ Jt n=-co^ ^ 
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(4.20) 


1, if j = 0 
0, if j ^ 0 


In actual computations the number of simultaneous equations 
must be truncated. One can, for example, solve a system of 
2N + 1 equations corresponding to -N ^ j ^ N. The choice of 
N must be such that the truncated version of the wave-number 
response, Eq. (4,3), do not change appreciably by further 
Increase of the number of terms used in the computation. 

Equations (4,20) are derived for finite and 

2 

finite K^. - Ico , and these equations cannot be reduced to 
those for supports rigid in translation or rotation. For 
example, if the supports are rigid in translation then Kj. 
becomes infinitely large, but the summation of all the 
must be zero since the deflection at each support is zero, 
and the product of these two becomes indefinite. To deal 
with this case, substitute 


Ao = - E A„ (4.21) 

^ n=-eo “ 
rv^O 

into Eq, (4,3) to obtain 


H(x,h) = Z A (e 
rr^O 


• Q 


) 


(4.22) 


Correspondingly, virtual displacements are chosen in the 
form of 


Skj (. 


) e 
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Then, instead of Eq,. (4.20), one obtains from a similar 
derivation 


J n*-oo« X X n=-«o X 

n¥0 n^O 


An 


« -1, j # 0 


(4.23) 


It is Interesting to note that if the supports are rigid in 
translation but without constraints in rotation (the case of 
hinge supports) the equations for Aj can be decoupled. For 
such a case, Eq. (4.23) reduces to 


AjCp(j) + cp(0)^|_^Ajj » -1, 3 0 (4.24) 

n^O 

Eq* (4.24) shows that the product AjCp(3) is independent of 3; 

1 • e • , 

A^Cp(l ) = A2Cp(2) = ..... (4,25) 

Thus , substituting 

Ajj = AjCp(3)/cp(n) (4.26) 

into Eq, (4.24), one obtains equations involving only one 
unknown! 




(4.27) 


which is solved readily to give 
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*3 - (*.28) 

Uncoupled solutions such as Eq.. (4.28) are not 

restricted to the case of hinge supports. In fact, if one 

of the two infinite stuns in Eq. (4.20) can be dropped, a 

substitution of the type of Eq. (4.26) is possible which is 

the key for reducing Eq. (4.24) to Eq. (4.2?) containing only 

one Aj in each equation. This is the case when either 
2 2 

K^-Mu> *0orI^-Iu) = 0ji.e., when the elastic supports 
offer no translational constraint or no rotational constraint. 

If the rotational constraint is infinite, one again 
cannot obtain a reduced equation from Eq. (4.20) which is 
valid only for finite support constraints. To derive a reduced 
equation one must use the zero slope condition at the supports 

Z = 0 (4.29) 

n=-oo 

The remaining procedure is very similar to that leading to 
Eq. (4.24) 

Now, we are ready to compute the spectral density of 
the structiural response. Once Aj are determined by either 
one of Eqs. (4.20), (4.23) or (4.28) depending on the problem, 
they can be substituted into Eq. (4.3) to obtain the wave- 
number response function, H. Further substitution of this 
H function into Eq. (2.36) along with the statistical proper- 
ties of the turbulence pressure gives the cross-spectral 



density of the structiiral displacement response 


k, 3 Numerical Example 


To Illustrate the application of the present theory* 
the spectral densities of the acceleration response at a mid- 
span location (i«e.» » x^ ” ^/2) have been computed based 

on the following physiceJ. datai 
properties of the beamt 

4 2 

D (bending rigidity) = 3.935 x 10 N-m 
i. (span- length) = 0,508 m 
m (mass per unit length) » 9.746 Kg/m 
properties of the. suzxounding fluid media i 
Pi ~ f*2 ~ P (density) = 0,11015 Kg/m^ 
a^ = a2 - a (speed of sound) = 26l,6 m/sec 
tJeo (free-stream velocity on upper side of beam) 

= 575.6 m/sec 

d (cavity depth) = 0,1178 m 
properties of the turbulent pressure C28]i 


_ , ^ 1 S 4 -K (|tu|S/Uoo) 

$ (0,o>) = spectral density = :5’ 2 A„ e « 

p ^ n=l 


tjf(S) = decay factor = exp(--^) 

o 

U„ (characteristic convection velocity of the 
c 

turbulence) = 0,75 Uoo 
S (boundary-layer thickness) = 0,279 m 


k6 

, , -2 
Aj^ = X 10 

^2 ^ 7.5 X lo"^ 

kj = -9.3 X lo"^ 

= -2.5 X 10'^ 

Except for the additional information about the span length X 
the above physical data are the same as those used previously 
for the unsupported beam, and they were taken from a recent 
experiment on a multi-panel system £28], The structural 
specimen used in this experiment was actually a two-dimensional 
panel array as shown in Fig. 4.2. Therefore, some data have 
been converted to their one -dimensional equivalents. For 
example, the bending rigidity 0 of the beam was the average 
value for the skin and the reinforcing stringers over a imlt 
width, and the specific mass m was obtained similarly. 

Although the actual structural specimen had only seven spans 
and the two end-spans were somewhat shorter, it was felt that 
the theory of an infinite periodic beam on evenly spaced 
supports should give a reasonable result for the response 
spectrum at the center of the middle span where accelerometer 
A20 was located (referred to Fig. 4.2), and where the effects 
of the end spans were least important. The translational 
constraints provided by the supporting frames were sufficiently 
strong to justify taking the translational spring constant ICj 
of the supports to be infinite (i.e., the deflections at the 
supports were assumed to be zero). For the rotational con- 

2 

straints we selected - 60 N-m/rad and I = 3*3 x 10 Kg-m , 
These are the one-dimensional equivalents of the torsional 


= 5.78 X 10 
Kg “ 2,43 X 10 
= 1.12 
= 11.57 




pig* ^*2 The experimental set-up 
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constraints of the frames if the torsional mode of each frame 
is a half-sine curve. 

In order to minimize the computer cost* the computation 
ivas carried out only to 3000 Hz and at the intervals of every 
50 Hz. The resolution of the computed spectrum was compromised 
somewhat hy the use of coarse intervals, but o\ir main objective 
was to find the general trend which could be revealed by the 
values at 50 Hz intervals. 

In Fig. 4.3 the computed spectrum is shown along with 
the experimental spectrum. As it is customary, the experimen- 
tal spectrum is one-sided (restricted to the positive frequen- 
cy domain) t therefore, the theoretical spectrum has been 
converted by multiplying the computed two-sided values by 
two. It also should be noted that the experimental results 
were obtained using a filter of 1 Hz bandwidth. This accounts 
for its much more rugged appearance than the theoretical one 
computed at much larger intervals of 50 Hz. Piarthermore, 
experimentally obtained signals may contain noise other than 
the structural response. Although the theoretical and the 
experimental curves show the same general trend, the former 
is lower than the latter throughout the entire frequency 
range investigated. This is to be expected since the theo- 
retical curve represents the average between the panel response 
and the stringer response, whereas the experimental curve 
shows the panel response alone. 



gVlte 


theoretical 



frequency I Hs 


■p- 


VO 


Pig* ^*3 Acceleration spectra of structural response 
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^.4 Conclusion 

The case of fluid- loaded infinite periodic beam* 
considered in this chapter, is one of the very few where 
a mathematically exact solution for the wave-niamber response 
function, H, can be obtained. If the beam is finite in length 
then the method of transfer matrix may be more preferable; 
however, the effect of fluid loading cannot be accounted for 
exactly (in the mathematical sense) at the present time. 

Further extensions to the two-dimensional case of 
panel systems are obvious. If only one row of panels is 
considered, and if the two parallel edges of the panel row 
are asstuned to be simply supported, then separation of spatial 
variables is possible in expressing the structural response. 
This is the well-known levy's type solution for plate 
problems. With small modifications, the solution for the 
one-dimensional beam case can be changed to suit such a panel 
row problem. When more than one row of panels are included 
in the structural model the separation of spatial variables 
in the structural motion is no longer mathematically exact, 
but a separable form can still be used as an approximation. 
Although new concepts are not required invtreating such two- 
dimensional problems, the machine computation time can become 
extremely excessive and burdensome to small research budgets. 
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V. MEMBRANE WITH RANDOM TENSION 

5.1 Introduction 

The unsupported infinite beam and the periodically 
supported infinite beam discussed in the preceding chapters 
are idealized models. However t actual structures cannot be 
constructed in ideal manners. The material properties vary 
randomly throughout the entire structure and manufacturing 
errors always exist in size and shape such as the span length 
and the cross-section of a beam, the pre-tension in a membrane, 
etc. It has been shown by Lin and Yang C 33- 35 3 that in the 
case of a periodic beam the randomness in the structure 
properties causes appreciable variations in the structural 
response from the ideal model results. In this chapter, the 
structural response of a membrane to the subsonic boundary- 
layer turbulence will be investigated. The pre-tension in 
the membrane will be treated as a random process in space, 
and the membrane is surrounded by an acoustic environment 
similar to that in the preceding chapters. 

5.2 Wave-Number Response Function of Membrane 

The structural, model chosen for the present study is 
a membrane which is infinitely long in the x-direction and is 
fixed along y = 0 and y = b as shown in Pig. 5.1. This 
membrane is backed on the lower side by a cavity of depth 
d which is filled with an initially quiescent fluid of density 
(*2 sdid sound speed 3 l 2 » ^ upper side the membrane is 



p(-turbulenoe prassxire) 



Figt 5«1 An infinite membrane with random tension under the 
excitation of boundary>layer turbulence 
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exposed to the excitation of a subsonic boundary-layer turbu- 
lent pressure p. The fluid on the upper side of the membrane 
which carries the turbulence has a free-stream velocity 
density and sound speed a^^. The design pre-tension in the 
membrane is Tq; however* it will be asstimed that this uniform 
tension can only be achieved in the x-directlon but in the 
manufacturing process a small rsuidom variation €f (x) has 
developed in the y-direction where € « 1, 

As the membrane responds to the excitation its motion 
will generate additional pressures in the fluid media on the 
upper and the lower sides. Denoting such pressures by Pj^ and 
P2t respectively* the governing equation of the membrane is 
given by 


m 


afw 

at^ 




- p + (Pi - P2)j=.a 


(5.1) 


where m is the mass of the membrane per unit area* Y is the 
viscous damping coefficient. 

As discussed in Chapter II* the turbulent pressure 
p should be replaced by expCi(tot - kx)] S(y - y*) and the 
structural response w by H(x*k*y,y* *u>) exp(iwt) for the 
purpose of determining the wave-number response ftaiction 
H(x*k*y*y' *co) , Thus* Eq, (5.1) becomes 

2 2 , 

C-mw^H + iroJH - Tq 

Sx Ssn 
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i(tut - kx) ^ ^ 

= e S(y - y) + (p - 


(5.2) 


The boundary conditions for the wave-number response function 
H are 


H(x,k,y,y*,u>) =0 at y = 0, b (5.3) 


These boundary conditions are automatically satisfied if H is 
expressed in a Levy type series 


H(x,k,y,y',a)) sin ( Aj^()i,k,y' ,to) d^ 

(5.4) 

In the sequel the symbol Aj^(^) will be used in lieu of Ajj(^, 
k,y',co) for compactness. 

For the determination of the radiation pressure 
we use the usual approximation that it can be calculated 
without regard to the presence of the txirbulence. Under this 
condition, the radiated pressure Pj^ is governed by 


2 2 2 

a 3 2 2 3 3 3 

+ W ”1 - ^ <7/ * ^ I'l “ ® ‘5.5) 


with the boundary condition 


3z 




ixijt 


(5.6) 


The exact solution for p^^ is unknown but in close proximity 
to the membrane p^ can be approximated by the expression 



55 


iojt - . ntiy 7 _jux ^ 

® nil J ® ^ <5.7) 

— oo 

Substituting Eq. (5.7) into Eq, (5.5) » multiplying both sides 
of the equation with sin(n*Try/b) exp(i^'x), and integrating 
over -00 < X < oo and 0 < y < b, one obtains 

g^(z,^) - gj^Cz.yu) = 0 (5.8) 

where 

=/ + (2^)^ - (5.9) 


Equation (5.8) is solved, with the condition that Pj^ only 
propagates in the negative z-domain, to give 


g^(z,y4) = Cjj(/x) exp[kj^(/^)zD (5.10) 

Substituting Eqs, (5.^)» (5.7) and (5.10) into Eq.(5.6), and 
using the orthogonality properties, we have 


®n‘/*> - Pi 


(<i> - 

lCn</^) 


^kn(/)d 


Ajj^) 


(5.11) 


Therefore, the radiation pressure p^^, when evaluated at 
z - -d, is 


z=-d 


= r, 


r 

n=l 


sin 


rory 


J 

— 60 


(uj - Uo^j^) 


*n</> 


(5.12) 
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The transmitted pressure is governed by the 
Helmholtz equation 


+ >L + »!, 

at^ ax^ 3y^ az^ 


) 


?2 =0 


(5.13) 


and is subject to the boundary conditions 

ap2 

— =0 at z = 0 (5.1^) 

2 iurt 

H e at z = -d (5.15) 

= 0 at y = 0. b (5.16) 


The conditions (5.1^) snd (5.16) corresponds to the ideal 
case where the bottom and side walls of the cavity are 
acoustically hard. The solution of this system is given by 


z=-d 



i(l}t eo oo 

e Z Z cos 

n=l r=0 

n+r=odd 


rrry 

T~ 


^ n 

7Tg(r) 


cot lCr.(/<)d 

kr(/) 




d^ 


(5.17) 


where 



(5.18) 

(5.19) 


Substituting Eqs, (5.^)» (5.12) and (5.17) into Eq. (5.2), we 
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have 


where 


-•o >o« 

- f~) S 94 - k) sin ^ 


- X)x 


dx3 dX 


( 5 . 20 ) 


v„(;^) " -«u? + iro) 4. T(,(^)* 4 To/i* - Pi 

<5.ai) 

. . X . 2 ^ 16 cot‘^(>A)d n i 

"1“ rSo yg(r) „2 - 

n+r^oda 

i+r»odd (5.22) 


In the present study Eq* (5.20) will he solved for the 
following two cases I 

case I t f(x) = 0 

M 

case II } f(x) = £ exp(ijux), Xq » 0 

Case I corresponds to the ideal uniform tension prohlemt i.e., 
the tension in the membrane is equal to the constant Tq in all 
directions. Case II represents the case where random varia- 
tion of tension in the membrane can be expressed as a super- 
position of sinusoidal variations. The fundamental wave 
number v will be assiuned to be detexministic and the complex 
amplitude St X^t random variables. 
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CASE I I f(xy = 0 


In this case Eq» (5*20) reduces to 

« 2 niry' „ 

V>*) * S(,/ - k) 

(5.33) 

which is an infinite set of simultaneous equations in Aj^. In 
actual computations the nximher of simultaneous equations must 
be truncated. One can solve a system of N equations corre- 
sponding to 1 ^ n £ N. The choice of N must be such that 
the truncated version of the structural response, (5»^)» is 
sufficiently accurate and it does not change appreciably by 
further increase of the niomber of terms. After truncating 
the number of simultaneous equations, Eq, (5*23) can be 
written in a matrix form 


CKi(jA)3 CA(|1)} = Cp} S(^ - k) (5.2if) 

where the elements of the matrix and the vector P are 

P = (h sin ^ (5.26) 

n D 0 


and is a Kronecker delta having the propertyi 


s 


nJ " 



n = J? 


(5.27) 


The solution of Eq, (5.2^) is 
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-1 


» C%(/a) 3 * £P} - k) 

Substitution of Eq, (5«28) into Eq., (5«^) yields 


(5.28) 


N 


nny 


H(x,k,y,y',tt>) » Z R sin 

n=l ^ ^ 


(5.29) 


vdiere Rj^ is the n-th element of the vector R which is defined 
as 

£r} = CKi(k)D“^ £P} e"^^ (5.30) 


M 


CASE II t f(x) '= Z Xi exp(ijvx)( 

3=-M *’ 


Xft » 0 


In this case Eq, (5.20) reduces to 




+ 


N 

Z 




J.M 


ill • 
(^) sin 


nuy* 

b 


S(^ - k) 


(5.31) 


where the number of the simultaneous equations has been 

truncated to N, In a matrix form 

M 

Ck^(;u«)3 £a(;a)} + eCKg] ^2 ^ Xj £a(^i + 3v)3 


» £p] S(jx - k) 

where the elements of the matrix K£ are 
nTT 2 <- 

'^rd ■ ^ni 


(5.32) 


(5.33) 
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which show that the matrix is a diagonal matrix. To solve 
Sq.« ( 5*32) is expanded in a series, keeping in mind 

that € is a small quantity, 

£a(;^)3 = C°a(^) 3 + eC^A(;4)3 + €^£^a(;a)3 + • • • (5.3^) 

Substituting Eq. (5»3^) into Eq, (5* 32), and grouping terms 
of the first and second powers of €, one obtains 

CKi(>a)3 £°a(^)3 = £p 3 - k) (5.35) 

1 M 0 

CKi(|i)3 £^a(;a)3 = - iKzl ^ Xj£ A(|i + jv)3 (5.36) 

2 “ 1 

£k^;(/^)3 £ a(;j)3 = - [%: ^ ^Xj£ A(ji + jv)3 (5.37) 
Eq. (5.35) is solved to give 

£°a(^)3 = CKi(;4)3“^ £p3 Kja - k) (5.38) 

Substituting Eq. (5.38) into Eq. (5.36), we have 




E X.£Ki(u + jv)3" 
js-M J ■' 


£p3 8(j^ + jv - k) 


(5.39) 


Likewise , 


£'"a(^)3 » CKi(/^)3"^CK23 ^e ^ ^^e^ x^ Xq^ 


M 


M 


£Ki(^ - qv)3"^[K23 CK^(/i + jV - qv)3"^ 



6l 


Cp 3 S(jA + jv - qv - k) 

2 

Thus* to the order of € , we have 

^ nTiy ** 3 

H(x,k,y,y*,aj) « Z sln-^CRjj + € Z ;F^ 

n=l ° 3=-M J " 

2 " “ * 3q 

+ € Z Z X. Xq G ] (5.^1) 

3=-M q=-M J ^ ^ 

where has been defined in Eq, (5.30), and and are 
the n-th elements of the following two vectors, respectively, 

Cp^ 3 = -CKi(k - 3 v)]"^CK 23 CKi(k)]"^£p} e"^^^ “ 

(5.42) 

CG^^3 = CK^(k - jv + qv)3“^CK2D [K^(k - 3v)r^CK2D 

C%(k)D’^P3 e-i^^ - (5.43) 

Now, if the random variables Xj, j =1,2,* * *,M, 
have zero means; i.e., 

ECX-3 =0, 3 = -M, • • •, -1, 1, • • *, M (5.44) 

3 

then the mean value of the wave-number response function H is 

eChD = z sin ^ £r^ + z z eCx. X ]gJ^3 
n=l ^ 3=-M q=-M •» q " 

(5.45) 

We note that the first term on the right side of Eq. (5.45) 
is the same as the wave-number response function for the ideal 
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uniform tension membrane « (5*29) » and the second term gives 
the contribution from the random variation in the tension 
field. The cross-spectrum of the response function H is 

ECH(x^,k,y^ H*(:^,k,y2,y^,(o)II 

N N nny i^ry ^ Vl Jt 

= E Z sin sin -tT-£R Ra + € E E 

n=l i=l ^ t> n X q=_jj 


EtXj CP^ p/ + R„ of + 4'* (5.46) 


where and y^ should be substituted into the variables with 
the subscript n and X2 and y^ should be substituted into those 
with the subscript Jt, 

For the special case where Xj are statistically in- 
dependent of each other and identically distributed random 
variables, and 


ECXj 



3 * q 
3 ^ q 


(5.^7) 


Eq,s. (5.^5) and (5.^6) reduce to, respectively. 


N 


ECh] = S sin ^ £r_ + -r— ^ .S 

n=l b 2M j=-M ^ 


y2 m2 M 


G^j} 


(5.^8) 


N N nny. Xiry2 ^ 

ECH H*] = 2: ^E sin sin — £Rj^ R^ 


2 2 

V^Tq “ 


J J 


3^* Jj 


2M 


jLm Vi ‘^n <5.49) 



where 
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(5.50) 


^ ijznr 
~ To 

is the coefficient of variation of the random tension in the 
y-direction. We note that due to the assiamption of uncorre- 
lated Xjf Eqt, (5.^9) becomes linear in each harmonic component 
of random tension field in the sense that the solution can be 
obtained by solving for each harmonic tension field vsuriation 
separately and then superimposed to obtain the total correla- 
tion. 

Now, the cross- spectrum of the structural response 
can be obtained by substituting into Eq. (2,52), Generally 
the integration over k in Eq, (2,52) must be performed numer- 
ically, but the integrations over y£ and y^* can often be 
calculated analytically. To illustrate we shall assume that 
the decay factor of the laressure field in the y-direction, 

~ expressed as 

4>2^yi " ^ 2 ^ = exp(- |y^ - y^|/Q) (5.5D 

where Q is an experimentally determined turbulence scale. 
Substitution of Eqs, (5,49) and (5*51) into Eq. (2.52) 
results in 




= *p(0,0,o>) 


W^(k - 


) £ 

c 


N 

L 

1=1 


mry. 


sin 


N N 
£ £ 
S=1 t*l 




Vt^To^ M 

+ -T 2 


sin 


b 
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+ V.j<*i> )33 13 ^, dk (5.52) 

where U^^g, ^ns,y "^ns,j elements in the n-th 

row euid s-th column of the following matrices, respectively. 


CO(x)3 = (f) CKi(k)]'^ 


(5.53) 


CVj(x)3 - - (-I-) CKi(k - Jv)3'^Cl%lCKt(k)] 
-i(k - jv)x 


,-l 


e 


(5.54) 


C»j(x)3 = (f) C%(k)3'‘CK2] C%(k - jv)3"‘CK23 

CKj(k)3-Ve-^'“ 


(5.55) 


and 


I 3 — ~S 

(l/Q)^ + (s7T/b)^ Q 

s7r/b tn/b 

* (1/Q)^ + (sVb)^ (1/Q)^ + (tVb)^ 

Cl - C(-D® + (-1)"^] + (-1)®'^‘^3 (5.56) 

In carrying out the numerical computations the most time- 
consuming part is the inversion of matrix beside the 
numerical integration over k which cannot be avoided. 
Therefore, in the actual computation, only the diagonal 
terms of Kj^ matrix will be kept to save the computer time. 
Hence the coupling of modes due to the effect of transmitted 
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pressure, pg, will be neglected, 
5.3 Numerical Example 


The following physical, data are used in the niomerical 
calculation! 

properties of the membrane (Mylar membrane)* 
f (density) » 1.39 x 10^ Kg/m^ 
t (thickness) = 8.89 x IG"® m 
m (mass per unit area) = (»t 
To (tension) = 61.29 N/m 
b (width) = 0.127 m 

Vj (coefficient of variation of random tension) = 0,1 
V (fiindamental wave number of random tension harmonics) 
s IT 


properties of surrounding fluid media* 

- Pz - density = 1.23 Kg/m^ 

Q-l ^ &2 ~ speed of sound = 340.36 m/sec 
l^(free stream velocity on the upper side of the 


membrance) = 30.48 m/sec 

d (cavity depth) = 2.54 x 10 ^ m 


properties of turbulence* 

- 1 ^ 

(0,0 ,u}) = spectral density = "o 7f~ E 

P ^ Uoo 



U (characteristic convection velocity of the 
o 

t\irbulence) = 0,8 Ueo - 1^1 

U 0 

i|fj^(S) = decay factor in thfe x-direction = e 
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y (IC - iiL) = — — ' 

1 Uc Tru^ec(l/Uj.e)^ + (k - (U/U^j)^l 

- JZLL 

= decay factor in the y-direction = e 
(Q = o<S*) 

S* (boundary- layer displacement thickness) 

= 0.001168 X Re“®*^ 

Re = V^x/v 
X = 0.609 m 

V = 1.464 X l0~^ mVsec 
and experimentary determined constants 
= 0.240 Kj = 0.470 

Ag = 1.08 Kg = 3.0 

= 1.80 = 14.0 

o( = 2 

6 (eddy lifetime) = -(1.24 x 10“^) (Uoo/a^) 

-3 

+ 1.15 X 10 sec 


Fig. 5.2 shows the spectral densities of structural 
response at Xj^ = Xg (the spectral densities of structural 
response is homogeneous in the x-direction, hence it is not 
a fxmction of x^^ or Xg) and “ y ^2 “ with or without the 
random variation in the membrane tension computed under the 
assumption of a frozen-pattern tiurbulence (i.e.» i|fj^(§) =1). 
Furthermore, only two harmonic terms of the random tension 
field, corresponding toV=7rm'*^ (j=l) and -V = -TT m ^ 

(j = -1), have been included in the computation. We recall 
that the cross-correlation of the structural response is 
linear in each harmonic component of the random tension field 
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when the random amplitudes of the tension* Xj* stre statisti- 
cally uncorrelated with each other. The spectral density of 
the structural response In the absence of random tension 
shows no peak within the frequency range shown In the figure, 
while several peaks appear If the random tension Is taken Into 
account. This change In appearance can be explained as 
follows. Since the exciting pres s\ure is assximed to be a 
frosen-pattern turbulence, the wave number k of the turbulence 
Is related to the frequency by k » Without random 

tension the wave number of the response, yn, of an Infinite 
membrane unsupported In the x-dlrectlon must be equal to 
Therefore, the coincidence resonance [3^3 would occur when one 
of the diagonal elements of the Influence coefficient matrix 
Ki co\ild become zero. If the Induced pressures, p^^ and p^, 
are neglected, this would happen when 


- mo) + T, 


(— )^ 
)'b ' 




'^O'U 


However, numerical computations have shown that the left-hand 
side of this equation Is always positive for any value of . 
Thus, coincidence can never occur when uniform tension Tq Is 
acting on the membrane alone. When the random tension terms 
are Included In the computation, there appear pertxirbatlon 
terms which have factors of the form of CK^(k - j = 

In the wave-number response function H. Thus 
a shift of the wave n\imber In the structural response from 
k to k - jv occurs. When j « 1, the coincidence frequency 
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can be calcxilated from the equation 

- nuo^ + - v)^ » 0 

c 

For n s 1 and 2, this equation has solutions 

f = 371 Hz, 720 Hz, for n » 1 

f = 425 Hz, 665 Hz, for n = 2 

No solution exists if n is greater than two. There is no 
coincidence resonance when j » -1, Note that these values 
are only rough estimates of the peak frequency since the 
induced pressure fields, p. and p , which have been ignored 
in the estimates provide additional inertia or stiffness to 
change the structural response, thus altering these peak 
frequencies. 

Fig* 5*3 shows the spectral densities of structural 
response using the measured spectrum of turbulence pressxire 
in the computation. Other specifications are the same as 
those used in obtaining Pig. 5*2. In this case, however, 
greater frequency intervals at every 50 Hz were used to save 
the computer time instead of 10 Hz interval in Fig. 5*2. 
Although the details are missing in this figure as compared 
to Pig, 5,2, the major effect of the random tension, which 
raises the spectrxm values at some frequencies by several 
orders of magnitude, can be seen clearly. The comparison 
with Fig, 5*2 also shows that the frozen-pattern assianption 
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is unconservative t resulting in gross underestimates in the 
structurail response calculation. 

5.^ Conclusion 

The ttirbulence decomposition scheme was applied to the 
two-dimensional problem of a membrane with random tension 
field. The Levy's series representation was used to describe 
the variation of the wave-number response function, H, in the 
y-direction. Therefore, integrations in the y-direction 
needed to calculate the cross-spectrum of the structural 
response could be carried out simply with only a single inte- 
gration remaining to be done numerically on a computer. The 
effect of the randomness in the structural properties was 
also investigated, in particular the random non-tmiform tension 
in the membrane. It has been shown that with a coefficient 
of variation of only 10^ in the tension field the structural 
response spectrum may be increased by several orders of 
magnitude at some frequencies. The variation of the tension 
field need not be accidental, but it may be caused intention- 
ally to create a spectral peak, for example, beneficial to 
skin-friction reduction. This proposition may be investi- 
gated in the future. The suitability of the frozen-pattern 
assumption was, again, examined in this chapter and the same 
conclusion as in Chapter III was drawnt namely, it leads to 
underestimation of struetursil response and omission of some 
important peaks in the spectral density of the response. Thus, 
the use of the frozen-pattern assumption shoiild be avoided for 
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the anal^is of the structural response spectra imder the 
boundary-layer turbulence excitation. 

In Chapter the emphasis has been placed on 

structural motions. In the next chapter our focus will be 
shifted to the skin- friction drag reduction itself. 
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VI. REYNOLDS STRESS ON A COMPLIANT SURFACE 
6.1 Introduction 

From observing the swimming of dolphins and amazed by 
the smooth motion of the fish Kramer hypothesized that 
favorable interactions between the flabby skin of the fish 
and water could reduce the skin- friction drag on the fish. 

Since then many researches, experimental [36-37!] well as 
theoretical [38-41 D, have been carried out in this field to 
find the mechanism of this highly applicable phenomenon. 
Aircraft and ship designers are especially Interested in this 
problem since reduction of skin-friction will result in a 
decrease in fuel consumption. It has been estimated that the 
skin-friction drag of an aircraft can be as high as $0f6 of the 
total drag. Although some experimental studies performed on 
very flexible membranes have indicated friction-drag reduction 
in fully turbulent boundary layers, the theoretical follow-up 
has not been as successful. For historical reviews of previous 
works the reader is referred to [36, 373* 

One requirement for a successful theoretical analysis 
is a thorough understanding of the complicated structural 
motion under the excitation of the boundary- layer turbulence 
[37] » to which Chapters III-V of this thesis were directed. 
Another requirement is the knowledge of the changing flow 
field resulting from the structural motion. In this chapter 
a perturbation approach similar to that used by Ffowcs 
Williams [482 snd Slick [4l3 will be applied to compute the 


7^ 


perturbation Reynolds stress near the structure-fluid interface. 
The analysis will be related directly to the two structural 
models considered in previous chapters. The first model, a 
one-dimensional unsupported beam will be used again to present 
the basic concepts, and then these concepts will be applied to 
the more realistic model of a one-dimensional infinite beam on 
evenly spaced supports. In both cases the turbulence decompo- 
sition scheme will be utilized. 

6.2 Theory 


Fig. 6.1 shows a one -dimensional infinite beam which 
is basically the same beam considered in Chapter III. As the 
beam responds to excitations its motion will generate addi- 
tional pressures in the fluid media on the upper and lower 
sides. As before, denoting these induced pressures by p^ and 
P 2 » respectively, the governing equation of the beam motion 
is given by 


m + Tj ^ + El jr 

at2 ' at 


A .2 

a w _ a w ^ , 

- T — o = P + (Pi - P 


ax‘ 


2^z=0 


Equation (6.1) differs from Eq. (3.1) in that a pre-tension T 
and a viscous damping coefficient have been included. 

The turbulence pressure field p can be expressed in 
the form previously introduced in Chapter II i 

ff i(ukt - kx) 


p(x,t) = 


dF(k.u) dG(u) 


( 6 . 2 ) 



p (turbulent pressure) 



Fig* 6,1 One-dimensional infinite beam under the 
excitation of boundary-layer turbulence 
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The meaning of the functions F and G has been discussed in 
Chapter II. We recall that in the integrand the circular 
i^equency co is related to the wave-number k and the component 
convection velocity u by the relationship u> = ku. Now, the 
induced quantities w, p^ and ^2 expressed as follows t 

00 

w(x,t) = II H(k.u) dP(k,u) dG(u) (6.3) 

— Oo 

00 

P^(x»y»'t) =11 Pj^(k,u,y) " ^^dP(k,u) dG(u) 

(6.4) 

Pg(x,y,t) = II P 2 (k,u,y) “ ^^dP(k,u) dG(u) 

(6.5) 

where H, Pj^, and may be called the wave-ntunber response 
functions for the structtiral displacement, the radiated 
presstire, and the transmitted sound, respectively. 

To determine the radiated pressure on the upper side 
of the beam, p^, we use the usual asstimption that the addi- 
tional pressure can be obtained without regard to the presence 
of the turbulent pressure p. Then, Pj^ can be obtained in the 
same way as in Chapter Illf 

2 

P. (k,u,0) = Pi^i i H(k,u) (6.6) 

where denotes the free-stream velocity of the flow on the 
upper side of the beam. For the case of incompressible fluid, 
we let aj^ -»"oo. Then Eq, (6,6) reduces to 
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Pj(k,u,0) 



H(k,u) 


( 6 . 7 ) 


Following the same procedure as in Chapter III» the pressxire 
generated on the lower side of the beam is found to be 


PgOc.u.C) 



cotC/u^ - a2^ kd/a2) 

- ag^ k/ag 


H(k,u) 


( 6 . 8 ) 


where the bottom wall of the cavity is assumed to be acoustic 
cally hard. Taking the limit ag -*■<>«>» we have the incompress- 
ible flow solution 


PgCk.u.O) = -fgU)2 ooth(kd) 

k 

For what follows we shall restrict our attention to the in- 
compressible case. 

Substituting Eqs. (6.2)-(6.5) ♦ . (6.7) f and (6,9) into 
£q. (6.1). one obtains 

H(k,u) = C-mu)^ + iTju) + Dk^ + Tk^ 

. ooth(kd) j~l 

This equation is similar to those derived in 

however, in these references the effects of the radiated 

pressure p^ was not taken into account. 

Let the velocity components of the turbulent boundary- 
layer flow on the upper side of the compliant plate be H t u' 
and V + V* . where "u and v" are the velocity components when 


pik(u - Uoo)^ 

(6.10) 
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the hovindary is rigid and u* and v* are perturbation velocities 
induced by the structural motion. Now assume that the pertur- 
bation velocities* u* and v'» can be expressed as power series 
in yi 

u*(x,y,t) = L fiL (x,t) y^ (6.11) 

n=0 “ 

v'(x,y.t) = Z b_(x.t) y^ (6.12) 

n=0 


From the conservation of mass of an incompressible fluid, 
we have 


au» 

dx 



(6.13) 


In order that this condition is satisfied the coefficients 
in Eqs. (6.11) and (6.12) must be related as followst 


b 

n 


(x,t) 


1 

n dx 


for n = 1,2,» • 


(6.14) 


Therefore, substituting Eq. (6.14) into Eq. (6.12), we have 


OO . (x,t) — 

v’(x,y,t) = bQ(x,t) - E — y*' 

^ n=l n ax 


(6.15) 


The boundary conditions at the beam sturface are 


u(w) + u'(x,w,t) * 0 


( 6 . 16 ) 


v(w) + v*(x,w,t) = 


3w(x,t) 

at 


(6.17) 


Taking Taylor expansions of u(w) and v(w) about y - 0, and 
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neglecting the high order terms of w, one obtains 


u(w) = u(0) + + . . . « (^)y^QW ( 6 . 18 ) 

v(w) = v(0) + (|j)y_QW + • • • « 0 (6.19) 

since u(0) and 'v(O) are zero and the continuity of the 

unperturbed fluid requires that 

- A 

'ay'y=o “ “ 'ax'y=o *’ ® 

Substituting Eqs. (6.11) and (6.18) into Eq. (6.16) * one 
obtains 

a^(x,t) + aj^(x,t)w + a 2 (x,t)w^ + • • • * ^ 

( 6 . 20 ) 


Keeping only the first term on the left-hand side of this 
equation, we have 

a^(x,t) = -Uw (6.21) 

where U is defined as 


U 


(lE) 

'3y'y=0 


( 6 . 22 ) 


In a similar way we obtain from Eqs, (6.15). (6,1?) and (6.19), 


l>Q(x,t) 


aao(x,t) 1 aa.(x,t) 5 

w - — — r w* - « • • « — 

ax 2 


(6.23) 
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Agaiiit retaining only the first term on the left-hand side, 

Equation (6.24) confirms what might have been concluded from 
intuitive reasoning that the perturbed velocity in the y- 
direction is equal to the structural motion in that direction. 
Substituting Eqs. (6.21) and (6.24) into Eqs (6.11) and (6.15) t 
we have 


u’(x,y,t) = - Uw(x,t) + S a (x,t)y^ 

n=l " 

3w(x,t) dw(x,t) 
V(x,y,t) = — + U y 


•f 1 n 

2 y 

n-2 h 3x 


(6.25) 


( 6 . 26 ) 


in addition to the continuity equation and the boundary 
conditions, the velocity components in the flow field must 
satisfy, two momenttua equations, one in the x-direction and 
the other in the y-direction. However, it is sometimes more 
convenient to replace one of these two equations by a diffusion 
equation gpverni^- the vortiqity perturbation! 


av' 

= aj- sy 


(6.27) 


In the proximity of the structure this diffusion equation 
is given by 


ato» 

at 


a^ a^ 


(6.28) 
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vdiere v represents the coefficient of viscosity. We shall 
apply this diffusion equation next to calculate the higher 
order coefficients in the power series of u* and v*. Since 
( 0 * results from the structiural motion it can he expressed 
in a form similar to Eqs. (6.2)-(6.5)t i.e.» 

oo 

f I Kuht - kx) 

u>*(x,y,t) « I n(k,u,y) e dP(k,u) dG(u) 

(6.29) 

Substituting this expression into Eq. (6.28), we have 

,2 2 iuk 

-^n(k,u,y) - (k + — )n(k,u,y) = 0 ( 6 . 30 ) 

Sy 

with the condition that the vorticity vanishes at infinity. 

A solution for this differential equation may be expressed as 


_ 2 i 

n(k,u,y) =Q(k,u) expC-(k + -^) y] (6.31) 


where the function Cl is to be determined. Expemding the 
exponential function on the right-hand side of Eq. (6.31) in 
a Taylor series about y=0, one obtains 


2 iuk A v*' n 

expC-(k + -v~)*yD = S -Jj- (-1) 

n»0 


(k" 


iuk 


_n 

.2 


Therefore, Eq. (6.29) can be written as 




OO' 

(x,y,t) =|Jn(k,u) + ^^)^ 


ei(ukt - kx) 


(6.32) 



However* from Eqs, (6,25)-(6,27) 


a^w(x,t) ^ ^ ^d^w(x,t) _ 

u)»U.y»t) = - a^(x,t) + C-^-^ — U 


- 2a,(x,t)]y + • • • 


(6.33) 


In consistence with the representations of w and oo'* Eqs. 
(6.3) and (6.29)* each a^ also can be expressed as follows t 


oo 

aj^(x*t) = j| Aj^(k*u) e^ 


(iikt - kx) 


dP(k*u) dG(u) (6.34) 


Substituting this equation and Eq. (6.3) into Eq. (6.33)* 


u)*(x*y*t) = 1 1 [uk^H(k*u) - Aj^(k*u) - Uk^H(k,u)y 


- 2A2(k,u)y3 dF(k*u) dG(u) 


+ higher order terms in y 


(6.35) 


n 


By equating coefficients of y in Eqs. (6.32) and (6.35) » 
we obtain 


n(k»u) = uk H(k*u) - Aj^(k*u) 


( 6 . 36 ) 


— 2 ixik i 2 

Xl(k*u) (k + -^p) * Uk H(k,u) + 2Ag(k*u) (6.37) 


Eliminating n from these two equations* we have 
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+ 2A2(k»u)3 


(6.38) 


To calculate and we use the momenttm equation 
In the x-dlrection 


3(u + u') a(u + u') a(u + u') 

- + (u + u’) + (v + V') 


at 


ax 


ay 


1 a(p + p. ) 2 

-- — — i- + v\r(u + U') 

Pi a* 


(6.39) 


However* for the rigid wall case, we have 


au ^ _ au . _ au i ap . , 2_ 

•~ + U-— + V = --i ^ + A)VU 

at ax ay Pi ax 


(6.40) 


Thus, by subtracting Eq. (6.40) from Eq. (6.39) and neglecting 
the square terms of perturbation velocities, one obtains 


au' 

at 


+ u 


au* 

ax 


+ u* 


!« + V ^ 

ax ay 


1 2 
+ V* u* 


au 

ay 


Pi ax. 


(6.41) 


It has been shown in Eqs^ (6.18) and (6.19) that at the 


proximity of the beam u and v can be ap^oximated as follows t 


u(y.) = u(b) ^ay^yaO y ■*' * * * “ 
av 

v(y) = v(o) + (^)y-o y + • • • = 0 

Thus, the Navier-Stokes equation for the perturbation 
velocities near the compli^t surface can be approximated as 
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au* 

at 


au* 


1 _ ^ 

Pi 


2 

W u* 


(6.42) 


Substituting Eqs. (6.4), (6.25), (6.26), and (6.3^) in Eq. 
(6Ji2) and letting y * 0, we obtain 

ikPj(k,u,0) = -v[Uk^H(k,u) + 2A2(k,u)3 (6.43) 

Therefore 

ik TT 2 

A2(k,u) = “ ---kH(k,u) (6.44) 

Eliminating A2(k,u) from Eqs. (6.38) and (6,44), 

A^(k,u) = uk^H(k,u) (k^ + ^)‘^P^(k,u,0) (6.45) 


Now, substituting Eqs, (6.3) » (6,44), and (6,45) into Eqs. 
(6,25) and (6.26), we have 


(( i(ukt - kx) 

u'(x,y,t) = 1 B(k,u,y) e dP(k,u) dG(u) 


— oo 


(6.46) 


v'(x,y,t) 


oo 

= j| C(k,u 


i(ukt - kx) 

,y) e dF(k,u) dG(u) 


(6.47) 


where 


B(k,u,y) * -UH(k,u) + £uk^H(k,u) + (k^ 




Pj^(k,u,0)}y + higher order terms in y (6.48) 
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C(kfU,y) » iukH(lCtU} - ikUH(lE»u}y 

+ higher order terms in y (6,49) 


The Reynolds stress associated with the perttirbation 
TelocitieSf u* and v*» is defined as^ 


•cCxty) = fiE£u*v*3 


(6.50) 


The us\ial minus sign associated with the definition of the 
Reynolds stress is dropped here since t as shown in Fig. 6.1 • 
the positive v* is taken to he opposite to the conventional 
direction. Now, substituting Eqs. (6.46) through (6.49) into 
Eq. (6.50) and applying £q. (2.26), one obtains 


'c(y) 


00 

= jj T(k,u, 


y) Sp(k,u) dk du 


(6.51) 


in which 


T(k,u,y) = /Oj^iukUHH* + fj^C-ikU^HH* - iu^k^HH^ 

^ uk^ ,2 . iukv-i„ 

+ hi^er order terms in y 


(6.52) 


Since co = ku zmd the wave -number spectriam of the turbulence 


# In theory, the imaginary i>arts in the expressions for u* 
8Uid V' should be zero* however, some approximations have 
been used in the analysis which may lead to non- zero 
imaginary parts in these expressions. Such superfluous 
terms should be discarded in calculating the Reynolds 
stress. Therefore, instead of Eq. (6.50) one can write 

T(x,y) a fj^ECRe(u') Re(v*)3 



86 


pressure t Sp, can be expressed in teznas of the frequency 
spectriun using £q. (2,33), SQ* (6,51) can be changed also to 


t:(y) *= j| T(k»Y»y) ®p(0fW) W(k - dk doo (6.53) 

— oo 

where $p(0,u)) denotes the frequency spectriun of the turbulence 
pressure and ? accounts for the decay property with respect to 
spatial separation. Due to statistical stationality and homo- 
geneity of the structural motion, the first three terms in the 
right-hand side of Eq, (6,53) vanish. Thus the expression for 
the Reynolds stress reduces to 


•C(y) = 


(uk 


^ (k^ + ^)“^Pi(k,^,0) H*(k,^).y 


$p(0,w) YOc - ^) dk do) 
c 


(6.54) 


We now define a Reynolds number, R, based on the 
wave convection velocity uj/k and the wavelength 27T/k; 


R 


(uj/k)(27T/k) 

2TTV 



(6.55) 


Assuming this Reynolds number, R, is large, the square-root 
in the integrand of Sq. (6.54) is expanded in a Taylor series 
about 1/R « Oj 


OC^ H. =. 



+ — r* 

* iR^ 



(6.56) 


Por values of R much greater than unity, the first term of the 
series alone is adequate. Thus, 
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ucik 

/ivuj 




$«(0,iu) Y(k - dk doo (6.57) 

c 

This expression of the Reynolds stress differs from Eq, (33) 
in C^83 in that the radiation pressure* Pj, takes place of the 
pressure fluctuation. As shown in the process of obtaining 
Eq. (6.4-2 )* the radiation pressure is more directly related 
to the perturbation velocities and, therefore, to the pertur- 
bation Reynolds stress than the unperturbed turbulence pres- 
sure on the rigid wall. Thus, this expression of the Reynolds 
stress seems to be more plausible than that in C^l,4-83, 

Now, substitutiong Eqs. (6.7) and (6.10) into Eq. 
(8.57) t and keeping only the real part which represents the 
physical Reynolds stress, we have 



« 1 2 
/avkul 



v2,n 2 

Ueo) £C-m<0 


4 2 

+ Dk + Tk 



2 2 coth(kd) 2 22 

U») - — k 1 + ->1 } y 


®p(0,co) ?(k - -gL) dk du) (6.58) 

c 

It is clear from Eq, (6,54) that the perturbation Reynolds 
stress vanishes when there is no pressure radiation} i.e., 
P^(k,^,0) 0. This happens when the wave propagation speed 

of the structural motion, u/k, becomes equal to the free-stream 
velocity of the flow, I^o* Of course, in the absence of struc- 
tural motion the perturbation Reynolds stress also must be zero. 
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If the beam motion is favorable the Rejmolds stress may become 
negative. It has been suggested C^8^ that such a negative 
Reynolds stress in the vicinity of the surface may deprive 
the turbulence in the boundary layer with energy supply and, 
therefore, the turbulence level may decrease. As seen in Eq. 
(6.5B) the sign of the perturbation Reynolds stress can be 
altered by changing the structural properties. However, the 
relationship is subtle and extensive numerical studies would 
be required to reach any quantitative conclusions. 


6.3 Periodically Supported Beam 

The structural model used in the preceding section 
for a preliminary investigation of the change in the Reynolds 
stress in the fluid due to structural motion was an unsupported 
infinite beam. However the results obtained there can be 
extended easily to the case of periodically supported infinite 
beam. This model is shown in Fig. 6.2, and it resembles more 
realistically the construction of an airplane fuselage. 

The governing equation of the beam motion not directly 
on the support is given by 

a^w aw 

m + “^ — + D(1 + ig) 4 - T 2 

a-fc* ' at av« 


= p + (p^ 



y=0 


(6.59) 



p (turbulent preseure) 






r iv. V 

1 

Pj^ (radiated pressure) 
— 


: H 


A 

^ 1 

^ ^ * 

^ 7777^ 

p^ (transmitted pressure) 





4 <p2-«2> 

/// 77 /)/ yyy ~ / ///// y // y ///// / 


Fig* 6*2 One-dinenelonal periodic beam under the 
excitation of boundary-layer turbulence 
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where rj denotes the viscous damping coefficient* g is the 
loss factor of the beam material, D is the beam rigidity, and 
m, T, p and p are the same symbols used in the preceding 
section. For the purpose of determining the wave -number 
response function, H(x,h,u), the ttirbulent pressure p should 
be replaced by expCi(ukt - kx)] and the structural response 
w equated with H(x,k,u) exp(iukt). Other induced quantities 
should also be replaced by their wave-number response functions 
in the same manner. 

Since the supports give rise to multiple reflections 
of the propagation wave in the structure, the same expression 
for the wave-number response function H as that used in Chapter 
IV is suitable: i.e., 

H(x,k,u) = 2 exp(-i^jjX/£ ) (6.60) 

n=— ^ 


where 




= I'q + 2nn = ki + 2n7T 
The radiated pressure pj^ is governed by the equation 


,a 3 -2 

^at " “1 


2 a^ a" 

a. (— + — ) p = 0 

ax ay ^ 


( 6 . 62 ) 


and subject to the conditions that p^^ can propagate only in 
the region y> 0, and that 


^ay ^y=o 


= P,(iku ^ H e^"^" 


( 6 . 63 ) 
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The solution of this system^ when evaluated at y=0, is given 

2 

oo K - Uo.) 


Pi = (^)C(u^ - 


i(ukt - J^nxA) 


(6.64) 


where 


“» ° “^4’n 


(6.65) 


The solution for the case of incompressible flow can be 
obtained by taking the limit a^^-veo) i.e. 


Pi “ „5 „ Pin ' 


i(ukt - }^x/t) 


where 


Pm = Pi(^) H„ 


( 6 . 66 ) 


(6.67) 


The transmitted pressure p^ in the fluid medium 
-d s y -c 0 is governed by 






( 6 . 68 ) 


and subject to the conditions 


>P2 

dy 


= 0 at y * -d 


(6.69) 


^P2 - , . iukt 

= -P 2 («k) H e 


at y = 0 


( 6 . 70 ) 
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The solution for P 21 when evaluated at y s o» is given hy 


2 oe i(ulct — u xXO 

P2 = f2(uk) 


Tn 


(6^71) 


where 


2 


2 //*n \ 2 


( 6 . 72 ) 


Taking the limit a 2 -»> oo» Eq* (6.71) reduces to the incom- 
pressible flow solution 


. , , ,2 - cothO^n^A) i(uirt - yi^x/1) 

P2 « -f2(uk) e 


(6.73) 


At this point it is of interest to note that the H 

n 

are not all independent of each other due to the constraints 
at the periodic simple supports. At these supports the 
deflections are zero. Thus by substituting x = 0 into Eq. 
( 6 . 60 ), 

H(0,k,u) = S ° 

n=-«o 

or = - Z H (6.7^) 

n^O 


Hence, an alternative representation for H is 


H(x,k,u) 


n=-oo « 

n!^0 


(6.75) 
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which is also a better representation from a computational 
point of view. Substituting Eqs, {6,66) t (6.73) and (6,75) 
into Eq, {6^*59) r one obtains 

C,(n) - ,(0) 

n?tO 

^ 6 , 76 ) 


where 


Cp(n) = -m(uk) + i‘)|(ulc) + D(1 










(6.77) 


To determine the coefficients we apply the well-lmown 
virtual work principle. Specifically, we assume a virtual 
displacement 


Sw ® SHj 






(6,78) 


The virtual work done by the internal and external forces 
must smn up to zero. Due to the spatial periodicity of the 
structureuL response and the virtual displacement, it is only 
necessary to apply the virtual work principle to one periodic 
unit. Thus, multiplying Sw to both sides of Eq, (6.76) and 
integrating over 0 < x < i, we have 


cp(o) Hj + cp(0) Z Hjj = -1 for j # 0 
n?fcO 


(6*79) 
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Equation (6.79) shows that the product qp(j) K. is independent 

>1 

of j; i.e. , 

Cp(l)l^ = Cp(2)H2 = • • • = Cp(j)Hj = • * • ^ 

Thus t substituting 

= HjCp( j)/cp(n) (6,80) 

into Eq, (6,79)* we obtain an equation involving only one 
untoiown 

HjCp(j)cp(0)^E = -1 for j 0 (6.81) 

which is solved readly to give 


H. = -Ccp(3)cp(0) E for j 0 (6.82) 

0 n=-oo cpCn; 


Denote the velocity components of the turbulent 
boundary- layer flow on the upper side of the beam again by- 
u + u* and v + v* where u and v are the components correspond- 
ing to a rigid wall. In view of Eq, (6,60), the perturbation 
velocities, u* and v*, can be expressed as power series in y 
as follows I 


u* (x.y.t) 



oo 

.S 

3 *-00 




n 


^i(ukt - /ijx/i) 


v'(x,y,t) 


E E 

n=0 j SS-.O0 




n 


^i(ukt - ) 


(6.83) 

(6.84) 


From Eqs. (6.83) and (6,84) and the continuity equation for 
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an incompressible fluid; 


au* . ay* 

» + - = 0 

ax ay 


(6.85) 


one obtains 


OO OO 


A A-'*?”" 

~ n-1 „ i(ukt - n-1 

“ S E ny B ^ e ^ *J ny 

n=l j=-oo “J 


Multiplying both sides of the equation by exp(i^^x/0 and 
integrating over 0 •< x c j!, we have 


B, 


’nj 


^n-1,3 » n = 1,2, • • • (6.86) 


Therefore, by substituting Eq. (6,86) into Eq, (6»,84) , 


VCx.y.t) = [Boj + i <^) An-1, jy"] 

i(ukt - /jX/?) 


(6.87) 


The non-slip condition on the plate sxirface can be 
stated as 


u(vr) + u'(x,w,t) = 0 


(6.88) 


and 


v(w) + V’ (x,w,t) = 


aw(x,t) 

at 


(6.89) 



Talcing Taylor expansions of u(w) and v(w) about y 0, and 
neglecting the hi^ order terms of w» we have 


u(w) = (0)yao ^ 




v(w) = 0 (6,91) 

Substituting Eq.s, (6,83) and (6^.90) into Eq, (6,88) „ one 
obtains 




, 2 N lOilct - ii\x/X) 

(Aoj + AjjW + A 2 iw"^ + • • •) e 




I H. 

'3y^y=0 j=_<so 3 


Keeping only the first term on the left-hand side of this 
equation, we have 


^03 ' «3 


(6.92) 


where 


Similarly, from Eqs, (6,87) t (6,89) and (6,91 )» 


(6,93) 


Boj ” iukHj (6,9^) 

Now, substituting Eqs, (6,92)- and (6,9^) into Eqs, (6,83) 
and (6,87) » we obtain 
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u»(x»y.t) « 




(6.95) 


V(x,y,t) « E (iukH^ - UH^y + S, ^ . .y”) 

3*-*» J ^ J n*2 n n-l,j*^ 

^i(ukt - /i^x/X) 


( 6 . 96 ) 


The vorticity perturbation 


w'(%.v,t) = 53r - sp- 


(6.97) 


must satisfy the diffusion equation neeir the beamt 


2 2 

aw* ,a a^ , 

— = v(--^ + — 2 ) 10 * 

at ax* ajT 


( 6 . 98 ) 


In view of Eqs. (6.95) and (6.96). we; may assiuae that the 
vorticity perturbation has the form 


10 


(x.y.t) . 1 n, ei<»« 

j=.oo 3 


(6.99) 


Substituting Eq. (6.99) into Eq. ( 6 . 98 ), multiplying both 
sides of the equation by exp(i^jx/i), and integrating over 
0 < X -<J, one obtains 


- C(^)" - 0 


(6,100) 


Under the condition that the vorticity can diffuse only in 
the positive y-domain. the solution of Eq. (6.100) is foiuid 


exp£-C(^)^ + ^^3^: 


y} 

Expanding the exponential function in a Taylor series about 
y = 0, and substituing the result into Eq, (6.99)» we have 

,-U.y.t)=. 1 n, z 


to 


5^.00 "d n=0 ^'Ji 

i(ukt - |ijx/L) 


e 




However, from the definition of the vorticity, to' also can be 
expressed as follows, using Eqs, (6.95) aiid (6,96), 

o)*(x,y,t) = £ £uk(~)H. - Aj_j - C(^)^UHj + 2A2j3y 

J— *0O 

+ higher order terms in y) 

( 6 , 102 ) 

Equating coefficients of y^ in Eqs, (6.101) and (6.102), one 
obtains 

(6.103) 

iljC(^)^ + + 2A2j (6.104) 


Eliminating D . from these equations, we have 
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Now -the Navier-Stokes equation for the perturbation 
velocities in the vicinity of the beaun can be written as 


Uy 

at ^ sx 


+ v'U 


i. IZl 

Pi.®* 



(6,106') 


Substituting Eqs, (6,66) t (6,95) aiid (6,96) into Eq, (.6,106) 
and letting y 0, we obtain 


Z 




PlV[U(^)'^Hj + 2A2 ;j]} e 


i(ukt - 


= 0 


(6.107) 


Multiplying both sides of the equation by exp(i^jX/£) and 
integrating over 0 < x < 1,. one obtains 


■^2j " " 2pjV ^ ^Ij 2 'I ' “j 


(6.108) 


Substitution of Eq, (6,108) into Eq, (6,105) yields 




1 

PlV 


rUl 

JL 


iuk^ -i,yx 

V J ^ X 


C(V-)" + (-V^)P. 
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(6,109) 


Now» substituting Eqs, (6,108) and (6,109) into Eqs, (6,95) 
and ( 6 , 96 ) t we have 


u'(x,y,t) 

v»(x,y,t) 


3 =-00 3 

•• i(ukt - 

Z S 4 e 
3 »-ao 3 


(6,110) 


( 6 , 111 ) 
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where 


R, = -UHj 




+ higher order terms in y 


( 6 . 112 ) 


S 




iukHj 



UHjy + higher order terms in y 


( 6 . 113 ) 


We note that the velocity components » u* and v'l obtained 
here are induced by a frozen- pattern component, expCi(ukt - 

oo 

kx)j, of the turbulent pressure p; therefore, .£ R. exp( 

oo 3“"** ^ 

and ^ exp(-i^jV^) are the wave-number response 
fxmctions of the velocity components, u* and v', respectively. 
Thus, the total velocities induced by the turbulent pressure 
p are obtained by superposition as follows t 

Oo 

u'(x,y,t) = jj gi(ukt - /»jX/i) dG(u) 

( 6 . 114 ) 


— Oo 
Oo 


v*(x,y,t) » J| dF(k,u) dG(u) 

( 6 . 115 ) 


The perturbation Reynolds stress near the beam can 
now be computed from 

t:(x,y) = f^EQu’V] (6.116) 

Substituting Eqs. ( 6 . 112 )-( 6 . 115 ) into Eq. ( 6 . 116 ), one obtains 



I 
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*c(x,y) 



£ £ 
JIB-OO js^OO 


T„j(lt.u.y) 


^i2rr(j - n)x/Ji 


Sp(k,u) dk du (6,117) 

where 


Tnj0c.u,y) « sj (6,118) 

As in the precedixig section, the integration variable u in 
£q, (6,117) may be replaced by o>/k and the wave-number spectrtim 
of the turbulence pressure. Sp. replaced by the frequency 
spectrum to give 


oo 

r(x.y) = jj ^£ 




i 27 r(j - n)x/Ji 


» (O.cu) ?(k - iii.) dk do) 

^ Ut> 


(6,119) 


where the specific expression for T^^j. obtained by substituting 
Eqs, (6,112) and (6,113) into £q, (6,118) is given as followst 

’■njO'.py) • Pii“OTn«j + 


+ higher order terms in y (6,120) 

Since the structtiral motion is a stationary random process in 
time, the first term in the expression of does not con- 
tribute to the integral when substituted into £q, (6,119), In 
this periodic beam case, however, the spatial homogeneity of 


\lL 
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the structiirai motion is destroyed because of the supports. 
Consequently} the perturbation Reynolds stress near a periodic 
beam is a function of the spatial variable x, A reasonable 
measure of the net effect of the perturbation Reynolds stress 
appears to be the spatial average defined as follows i 

X 

<t:(y)> = J I ■c(x»y) dx (6,121) 

0 

Again, define a Reynolds number 

( 6 . 122 ) 

Eqs. (6.119) and (6.120) may be substituted into Eq. (6.121), 
and the result is simplified by expanding the square-root 
term in a Taylor series about 1 /r = 0. Keeping only the real 
part of the expression which represents the physical average 
Reynolds stress, one obtains 


<r(y)> = II /2v|u)| ^ ^ 


$ (0,0)) Y(k - ^) dk d(t) ( 6 . 123 ) 

P Oc 

Finally, substitution of Eqs, (6,6?) into Eq. (6.122) yields 


Z 

n=-oo 


$ (0,(o) W(k - ^) dk do) 
P ‘^c 


(6.12if) 
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6.4 Concluding Remarks 

Expressions for the perturbation Reynolds stress 
Induced by fluid-solid interaction were obtained for two 
structural models. The first structural model, an Infinite 
unsupported beam, was used to develope the basic concepts of 
the analysis. The results were then extended to the second 
model of an infinite beam, simply-supported at equal intervals. 
The second model is a more realistic model for the typical 
fuselage construction of an airplane. The effect of the radi- 
ation pressure was included in the formulation. The pertur- 
bation Reynolds stress involves a double integration with an 
integrand depending on the beam motion and radiated pressure. 
The complexity of the expression permits only qualitative 
discussions in this thesis. Hopefully, quantitative results 
can be documented after extensive numerical studies in the 
futtire. 

The radiation pressure, p^, required in the present 
analysis has been obtained from a wave equation for inviscid 
fluid. More rigorously the viscosity in the fluid and the 
velocity profile in the boundary layer should be taken into 
account when determining the radiation pressinre. However, at 
the present time no closed form solution is known of this more 
accurate wave equation. 


li 
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VII. GENERAL CONCLUSION 

This investigation has been concerned with the inter- 
action between a turbulent flow and certain types of structure 
responding to its excitation. The turbulence is typical of 
those associated with a boundary layer « having a cross-spectral 
density indicative of convection and statistical decay. It 
has been shown that a decaying turbulence can be constructed 
from superposing infinitely many components, each of which is 
convecting as a frozen-pattern at a different velocity. This 
turbulence decomposition scheme reduces greatly the computation 
time by reducing to one-half the number of integration which 
must be performed on a computer. Furthermore, the scheme 
provides a convenient way in which experimentally measxired 
cross-spectral density of the turbulent pressure fluctuation 
can be incorporated directly into the computation. 

The results of the structure-tiirbulence interaction 
were presented in terms of the spectral densities of the 
structural response and the perturbation Reynolds stress in 
the fluid at the vicinity of the interface. A number of 
structural models were considered in the investigation. 

Among the one-dimensional models were an vinsupported infinite 
beam and a periodically supported infinite beam. The first 
model was used to develope the basic ideas which were then 
applied to the more realistic second model resembling the 
fuselage construction of an aircraft. For the two-dimensional 
case a simple membrane was used to illustrate the type of 
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formulation applicable to most two-dimensional structvires. 
However ( a small random variation in the membrane tension 
was included in the analysis since ideally uniform tension 
never exists in practice. Moreover! the mathematical approach 
used in dealing with random membrane tension can be adapted 
to treat other random structural properties in general. Both 
the one -dimensional and two-dimensional structures mentioned 
above were backed by a space filled with an initially 
quiescent fluid to simulate the acoustic environment when 
the structure forms one side of a cabin of a sea- or air- 
craft. 

It has been found that important spectral peaks of 
the structiural response will not appear if decays in the 
tiarbulence is neglected in the analysis. Thusi the usual 
Taylor's hypothesis of frozen- pattern turbulence is uncon- 
servative as far as the assessment of structural reliability 
is concerned. The perturbation Reynolds stress is indicative 
of the change in the skin-friction drag due to structural 
motion. It has been shown that, given the statistical 
information of the boundary-layer turbulent pressure field, 
the perturbation Reynolds stress can be altered by varying 
the structural parameters. Therefore, the present study is 
potentially useful for designing flight or marine structures 
to minimize the total skin-friction drag. 
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